DEPARTMENT OF MATHEMATICS 


HASA OR b 6 8 3 

CLEMSON UKTI-VEIRSIT'Y 


OL-EJMSOIST, SOUTH OA-HbOLiIlsrA 


August 



Office of University Affairs, Code Y 
National Aeronautics and Space Administration 
Washington, D. C. 20546 

Gentlemen : 

Subject: Final Report, Grant NGR 41-001-016 

In accordance with NASA Research Grant NGR 41-001-016 
to Clemson University the following is the final report . 

Numbers enclosed in brackets refer to the List of References 
at the end of the report. 

The objective of the research sponsored by this NASA 
Grant was to study the existence, stability and computation 
of periodic and almost periodic solutions of vector differential 
equations x = Ax + F(x,t) and the analogous functional 
differential equations. The mathematics used for such systems 
treats this system as a perturbation of the linear system 
x = Ax. Those cases where the construction of almost periodic 
solutions leads to a trigonometric series with coefficients in 
which arbitrarily small divisors occur were of particular interest 
For convenience we refer to -these cases as the "small divisor" 
problem. 

The existence, stability and computation of an almost 
periodic solution is well known [4] if the perturbation F(t,x) 
is small enough and either the eigenvalues of A have nonzero 
real parts or the eigenvalues of a certain Jacobian matrix have 
nonzero real parts. Other cases were, generally intractable 
except by using the relatively crude methods which are known 
for the "small divisor" problem. The methods [1] for the small 
divisor problems are limited to very small perturbations and 
are necessarily troublesome to execute. 

We have discovered during this investigation a very 
significant generalization of the problem. It is possible in 
many cases to solve a nonlinear initial value problem 

x = f(t,x), x ( t ) = y 

for all y in some region of R n . It is then possible to study 
the existence, stability and computation of almost periodic 
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or periodic solutions of perturbed nonlinear systems 
x = f(t,x) + g(t,x). Many of the theorems ..concerning the 
perturbed linear system may be proved for the perturbed 
nonlinear system. One obvious benefit of such theorems is that 
previously the internal nonlinearities of such systems were 
treated as noise; and consequently small external perturbations 
would nullify the existence of the desired periodic or almost 
periodic solutions. Another benefit is that this technique 
provides a method to handle many problems previously tractable 
only by small divisor methods. This approach to perturbed 
nonlinear systems was discovered independently by L. E. May [5]; 
however, his analysis requires that the unperturbed system is 
essentially a linear system. Two papers illustrating the new 
approach to perturbed systems have been written and are described 
below, (4), (5). Work is progressing on other aspects of this 
technique, and new papers will be submitted soon. 

The work on "small divisor" problems was somewhat disap- 
pointing.. One paper was printed as a NASA contractor's report 
and is described below. In addition we have developed a formal 
recursive procedure for constructing almost periodic solutions 
in some situations; however, we have been unable to prove the 
convergence of the series generated by this procedure. The work 
which was done to establish the size of allowable perturbations 
was very complicated and it was found that in general the 
allowable perturbations were extremely small. 

The existence, stability and computation of periodic and 
almost periodic solutions to functional differential equations 
is analogous to that for the ordinary differential equation. 

One paper was submitted concerning the problem; however, much 
work is left undone. In particular, the generalization to 
perturbed nonlinear systems mentioned above for ordinary 
differential equations has not yet been done; however , no 
particular difficulties are forseen for this case. 

The following papers have been submitted for publication 
as a report of research performed under this grant. 

(1) "A family of solutions of certain nonautonomous 
differential equations by series of exponential functions" by 
T. G. Proctor and H. H. Suber. NASA Contractor Report 
NASA CR-1294 . This paper presents two theorems concerning the 
computation of periodic or almost periodic solutions of a 
differential equation. The first result makes use of a technique 
developed by Golomb [ 3 ] to construct a periodic solution to a 
nonlinear periodic differential equation, x = f(t,x), f odd in t, 
using a recursive technique to find the coefficients in a 
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trigonometric series. The technique is also applicable to 
quasiperiodic equations; however, we were able to guarantee 
convergence of the trigonometric series only in the periodic 
case. The coefficients can be calculated using a computer; 
however, we found the process time-consuming. The second 
result gives the existence and construction of an almost 
periodic solution to a differential equation of the form 
x = f(t,x) where f is odd and almost periodic in t and is small. 
This problem has the distinction of being the easiest small 
divisor problem in differential equations which has appeared in 
the literature. The paper is written so that the reader can 
determine the size of the allowable perturbation function 
(in this case f (t ,x) ) ; however, even in this relatively simple 
problem the resulting calculations are quite time-consuming. 

The two theorems described above are then applied to .the problem 
of reducing a 2-vector linear differential equation y = P(t)y 
to the easily handled form z = Az . 

(2) "Uniqueness and successive approximations for functional 
differential equations" by T, G. Proctor. This paper will appear 
in the Journal of Mathematical Analysis and Applications . The 
theory of functional differential equations is not as complete 

as that of ordinary differential equations. This paper was 
written to provide a basis for studying the existence and 
computation of solutions of functional differential equations 
possessing certain properties, for example, periodicity. The 
first result establishes a differential inequality analogous to 
the Kamke differential inequalities for ordinary differential 
equations. The use of such a theorem is illustrated repeatedly 
in (5) below. The theorem is then applied to the existence and 
construction of solutions to the initial value problem 

y(t) = <j>(t), a < t <_ t 0 , 

y(t ) = F(t ,y ( • ) ) , t Q < t < t Q + a 

where <p (t ) is a given function and F is a delay functional. In 
a later paper we plan to use the analysis above to study the 
existence and construction of periodic solutions to a periodic 
functional differential equation. 

(3) "Characteristic multipliers for some periodic 
differential equations" by T. G. Proctor. This paper will appear 
in the Proceedings of the American Mathematical Society . This 
paper gives two new theorems which give the construction of 
periodic solutions to nonlinear differential equations. The 
smallness condition on the perturbation function is replaced 

by finding an upper and lower trial solution. These theorems 
are quite useful in constructing examples of systems which 
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possess periodic solutions in systems where the perturbation is 
not small. The same construction is valid for almost periodic 
differential equations; however, the hypothesis is hard to 
verify in this case. The results are applied to the problem 
of finding characteristic multipliers for the 2-vector linear 
differential equation x - P(t)x, P with period T. A knowledge of 
the characteristic multipliers and vectors enables one to 
determine the solutions for all t by calculating the solution 
values for 0 <_ t £ T. 

(4) ’’Periodic solutions for perturbed nonlinear differential 
equations" by T. G . Proctor. This paper has been submitted to 
the Proceedings of the American Mathematical Society , Here we 
generalize the basic results in periodically perturbed linear 
systems to periodic perturbed nonlinear systems. The results 
should give new information concerning the motion of nonlinear 
springs, etc., in the presence of noise. The advantages of 

such a theory have been pointed out above. The chief difficulty 
comes in solving the unperturbed nonlinear initial value problem; 
however, this is known in many applications. The results can 
also be applied to the determination of characteristic multipliers 
of periodic 2-vector systems x - P(t)x as in (4) above. Here 
one reduces the problem to a Ricatti equation with periodic 
coefficients. If these coefficients are nearly constant the 
problem is easily solved since solutions to Ricatti equations 
with constant coefficients are well known. 

(5) "Integral manifolds for perturbed nonlinear systems" 
by H. H. Suber. This is Dr. Suber’ s dissertation and has been 
approved by Clemson University. Currently Dr. Suber is 
condensing the results and will submit them for publication soon. 

The main theorem in this paper generalizes basic results 
concerning the existence of integral manifolds for perturbed 
conditionally stable linear systems of ordinary differential 
equations [2] . Sufficient conditions are given for the exis- 
tence of an integral manifold for perturbed nonlinear system, 
in a neighborhood of a critical point, periodic orbit or periodic 
surface, in the case where the unperturbed state exhibits 
conditionally stable behavior. Although the proof of the 
theorem is constructive in nature it is not expected that the 
technique will be of great value for computational purposes. 
However, estimates for allowable size of perturbation are 
readily available. 

One copy of each of the papers above is enclosed as part 
of this report. 
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Dr. Suber's final two years were sponsored under this grant 
during which time he contributed papers (1) and (5) above. 

Talks covering the papers (1), (2) and (3) were presented to 
the Mathematical Association of America (Spring 1968) and the 
American Mathematical Society (Fall 1969 ). Three one-hour 
presentations concerning the papers above were given at Langley 
Research Center during the period of the grant. 

The work on this project may be summarized by saying that 
we have established several techniques for proving the existence 
and for computing periodic solutions for perturbed linear or 
nonlinear differential and functional differential equations. 

One small divisor problem has been treated] however, no real 
progress has been made on this class of problems. 


Prepared by 

T. G. Proctor 
Associate Professor 


cc: Dr. E. S. Armstrong (2) 

NASA Langley Research Center 

Mr. Mel ford A. Wilson 

Vice-President for Business and Finance 
Clems on University 

Enclosures: (9) copies 
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UNIQUENESS AND SUCCESSIVE APPROXIMATIONS FOR 
FUNCTIONAL DIFFERENTIAL EQUATIONS* 


T . G. Proctor 

Mathematics Department 
Clems on University 
Clemson, South Carolina 29631 


1 . INTRODUCTION 


Differential equations which express y'(t) as a function of 
past and present values of y (t ) have been called delay functional 
differential equations. Many of the ideas from the theory of 
ordinary differential equations have been generalized for this 
type of equation including the basic ideas of existence and 
uniqueness of solutions of initial value problems. 

Let a and t be numbers where -°° < a < 06 and t (> a) is a 
o — o — 

finite number. In case a = - 00 read [a, t] as (a, t ] . Suppose 
that 4> ( t ) is a prescribed continuous n-vector function on [a, t o ] 
and we wish to find a continuous function y(t) on some interval 
[a, t Q + a], a > 0, such that 


y (t ) = 4> (t ) , 
y’ (t) = F (t j y ( • ) ) , 


a <_ t <_ t Q , 

t < t < t + a, 

o — o 3 


( 1 ) 


where F(t, i|>( • ) ) is a function (functional) defined for t in 
[t Q , t + a] and ip in C(t) and taking values in R n and where 
C(t) is all continuous functions from [a, t] into some set D 


in R n . It is known [1] - [ 3 ] that for appropriate sets D if F is 
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continuous and satisfies a Lipschitz condition in ip, then the 
initial value problem has a unique solution and this solution 
may be constructed by successive approximations. 

In this paper we give theorems analogous to the Nagumo and 
Wintner theorems for ordinary differential equations [4] which 
prove the solution of the initial value problem (1) for other 
conditions on P is unique and can be constructed by successive 
approximations . 


2. A DIFFERENTIAL INEQUALITY 

First we give a lemma on differential inequalities which is 
analogous to the Kamke differential inequality theorems for 

ordinary differential equations [4], [5]. Let D be a region in 

n n 

R and let F be a function taking values in R for t <_ t <_ t + 

and 4 1 in C (t), where C^Ct) is all continuous functions from 

[a, t] into D, with the properties that lim F(t , i|> n ( • ) ) = 

F(t, ip( • ) ) whenever t, t^, •••are in .[t , t + a], ip, ip 1 ••• 
belong to ( t Q + a) and t t, 'ip ■* ip (in the sup norm 
topology). (If a = C]-) (t ) is all continuous functions from 

(- 00 , t ] into some compact subset of D.) Further , let (t, ip^) , 

(t, ipp ) in the domain of F and ip^is) <_ ^2 ^ s ) Ton a < s < t imply 
F(t, ip^( • )) £ F(t, \p 2 ( • ) ) . Here a vector inequality < means that 
every component of the left vector is less than the corresponding 
component of the right vector, etc. Also note that we use the 
same notation for a function \p given on an interval [u, t] and 
its restriction to a subinterval [a, t'] C [a, t]. 



3 


LEMMA 1. Let z be a continuous function from [a, fc + a] into D 
and be such that 

D + z(t) < F(t 3 z( •) ), t Q < t < t Q + a; (2) 

and let y ( t ) be a continuous function with the properties 

(a) y (t ) > z(t), a < t < t Q , 

(b ) y ' (t ) = F(t , y ( * ) ) , t Q £ t £ t Q + a. 

Then y(t) _> z (t ) for t < t £ t Q + a. (Here D + z(t) = 
lim sup[z(t + h) - z(t)]/h as h ->• 0 + .) 

Proof. Suppose y(t) £ z(t) on some largest interval [t o , 6] in 
[t Q , t + a]. (Here 8 could be t .) If for some 
k e {l, 2, •**, n) ,z k (6) = y k ( 5 ) then 

D + z k (6) < F k ( 6 , z ( • ) ) < F k ( 6 , y ( • ) ) = y k (6) , 

which implies z(t) <_ y ( t ) for values of t in a right neighborhood 

of 6, thus 6 = t + a. 

’ o 

o/i| 1/2 

The monotonic behavior of F is crucial since y(t) = 2 t 
is a solution of y ' (t ) = 1/y (t/2 ) on 1 £ t < and z(t) = t^ 

satisfies (1) on 1 £ t < 4 / ; however, z(t) > y(t) for 

2 3/ l 4 < t < . Thus, the Kamke differential inequality type 

theorem requires additional conditions. The continuity condition 
on F is sufficient to guarantee the existence of at least one 
solution to (1), see [1], [2] . The differential inequality theorem 
above then will infer the existence of a right maximal solution to 
(1) as Coppel [5] does for ordinary differential equations. Then 
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the strict inequality in (2) may be replaced with < provided y(t) 
is the maximal solution of (1), see Coppel [5] . 

3. A UNIQUENESS THEOREM. 

Let 4) be a given continuous n vector valued function on 

[a, t Q ] and let C^(t) be the collection of continuous functions 

f(s) on [a, t] which agree with <p on [a, t ] and which lie in 

S(4>(t o ) , b) for t Q < s < t < t Q + a, (S(y Qa b) = {y: | y-y Q | < b}). 

Let F be a function taking values in R n for t in [t , t + a] 

and ip in C b (t) with the property that lim F(t , ip (•)) = F(t, ip ( •)) 

" J n-*- 00 ^ ^ 

whenever t, t ±> •••, are in [t Q , t Q + a], ip, ip ± , •••, belong to 

C k (t 0 + a) and t n t, ip n -*• ip (in the sup norm topology). 

In addition to the continuity requirement above, we want to 
impose a restriction on F which is analogous to the Nagumo con- 
dition. For this purpose we define a set of linear functionals. 

Let W be the set of all linear functionals w, of the type 

t +a 

■ o 

w ( t , x ( • ) ) = x(s)d n(t, s), 

J a s 

defined for t Q < t <_ t Q + a, x a continuous function on 

[a, t Q + a] where ri(t, s) is as described below and where the 

functional differential equation 

x(t ) = 0 , a < t < t , 

° (3) 

x(t)' = w ( t , x ( • ) ) , t Q < t £ t + a, 

has a solution f (t) on some interval [a, t + <5 ] , where <5 >0, 

w o w w * 

such that f* w does not vanish on (t Q , t Q + 6 ] and 

li m ( ^ ) / 1 - t / 0. Here r)(t, s) is 

t+t + w 0 

o 
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on 

(i) defined., real-valued/(t 0 , t + a] x [ a, fc + a] 
(ii) constant for t £ s, 

( .ill ) nondecreasing in s, 

(iv) continuous in t uniformly with respect to s. 


( 4 ) 


J 


Some examples of such linear functionals w in which f (t) = t p , 

w 


0 < p £ 1 are 

(a) w(t, x ( • ) ) = 


t-e 


x(s )ds where a = -e, t = 0, 6 = e. 


(b) w(t, x(*)) = z-[Ax(t) + Bx(t/2)] where a = t = 0, A, B, 

U O 

such that p = A + B/2 P , and 
1 2 

(c) w(t, x(*)) = — p x(t ) where a = t = 0. 

t^ ° 


THEOREM 1. Let F be as specified above and satisfy 


I F ( t , - F(t, 

for (t, (t, ^ 2 ) in the domain of F, for t > t and for some 

w in W, Then the initial value problem (1) has at most one 
solution on any interval [a, t + e] for small e > 0. 

Proof, Assume there are two different solutions y ^ (t ) and y 2 ( t ) 
on some interval [a, t + e] . By shrinking e we may suppose that 
for y(t) = y-, (t ) - y Q (t ) , |y(t^ + e ) | 0 and e < 6 . Consider 

the solution x(t) = 0, a < t < t x(t) = i|y(t + e)|f ( t )/ 

| f w ( t o + e ) | , t < t £ t + e of the initial value problem (3). 

By (5) we have 

D + |y(t ) j £ w(t , |y ( • ) I ) , t Q < t £ t Q + e. 
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For any 0 < a < e the initial value problem 

{ 0, a < t < t , 

x(t ) , t Q <_ t <_ t Q + a , 

z ? (t) = w(t , z (t ) ) , t + a < t <_ t Q + £ , 

has a unique solution, namely x(t), since w(t, x( • ) ) satisfies a 

Lips chit z condition in x for t in [t + o, t + e] . The remarks 

following Lemma 1 then imply that we cannot have x(t ) > |y(t)| on 

a right neighborhood of t Q . Hence, there is a decreasing sequence 

of points t, t such that 
k o 

|y(t k ) | > x(t k ) , k - 1, 2, 

Thus “ t o ); ->• 0 but this is impossible since 

11m f w (t)/t - t Q ? 0 as t -*■ t o + . Hence y(t + e) = 0 . 

Remark 1. The proof of the theorem makes use of the linearity of 
the functional w and the form of the solutions of (3)> It is not 

readily apparent how one would generalize Kamke ’ s uniqueness 

theorem [4] since if w(t, x(»)) is a nonlinear functional, a 
solution of (1) passing through ^-|y(t + e) | may not exist. 

Remark 2, It is possible to relax the requirements on the 
integrator function n(t, s) in the case of nth order scalar 
functional differential equations as shown by Theorem 3 below. 
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4. SUCCESSIVE APPROXIMATIONS, 

Next we wish to show that the sequence of successive approxi- 
mations given by 

y n (t) = <Kt ) , « < t < t Ql n - 1, 2, ••• 3 

y x (t) = *(t 0 ) , (6) 

f t I ‘o ' 1 * f ’ 

y n (t) = 4> (t Q ) + F(s , y n _ 1 (* ) )ds , J n = 1, 2, • • • > 

t o 

converges for t in [a, t + 3] to a solution of (1) where 
t) 

3 = mlnja, 6^} and where M is such that | F( t , C • ) ) | £ M when 
t Q £ t £ t Q + min{a, <$ w } and \p is in C, (t + minja, 6 W }). It is 
easy to see that (6) defines a sequence, 

THEOREM 2 . Let F be as specified above and satisfy 

I F(t , ^ 1 (*)) - F (t, ^ 2 (*))| < w(t , I^-l - ^2 I ( * ) ) > 

for (t, ij^), (t j ^ 2 ) in the domain of F and t > t Q , for some w 
in W. Then the sequence (5) converges to the solution of (1) on 
Ct 0 , t 0 + 3], 

Proof, The sequence (6) is uniformly bounded and uniformly equi- 

continuous on [t , t + 3] so by Ascoli's theorem there is a 

o o 

subsequence {y (t ) }°° converging to some function y(t ) uniformly 
n k k=l 

on [t o , t Q + 3], Hence the sequence {y n + ^(t ) } converges 

uniformly to some function y*(t) on [t Q , t Q + 3] . If we can 
show A (t ) = 0, where 

A ( t ) = 11m sup ] co ( t ) | j w (t) = y (t) - y ,(t) 

n oo 11 11 11 • 


5 
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then y = y* = the unique solution of (1) and thus the whole 
sequence {y n (t ) J n=1 converges to y(t) uniformly on [t Q , t + 6] 

(see Hartman [4], pp. 4 and 4.1). The remainder of the proof 
amounts to first showing D + A(t) < w(t, A ( t ) ) . Then if we suppose 
A (t ) ? 0 we may repeat the strategy in Theorem 1 to get a con- 
tradiction. The details in the case of ordinary differential 
equations are given in [6] and [7] • 

5. A UNIQUENESS THEOREM FOR SCALAR PROBLEMS. 

Let $ be a real-valued function with n-1 continuous deri- 
vatives defined on [a, t ]: and for t in [t , t + a] let C, (t) 
be the set of all real-valued functions f with n-1 continuous 
derivatives defined on [a, t ] which agree with cj> on [a, t ] and 
satisfy 

|f (k) (t) - <J> (k) (t 0 )| < b, k = 0 , 1 , • • • , n-1, 

t in ( t Q , t + a] , (f^ k ^ is the kth derivative of f ), Let 
U(t, x ( • ) ) be real- valued for t in [t Q , t + a] and x in C b (t ) 
and be such that lim U(t , ^ (•)) = U(t , ( • ) ) whenever 

n-s-oo 

t, t la t 2 , •“ is in [t Q , t Q +a] , **• belong to C b (t Q + a) 

and t n t, \p n \p (in the sup norm topology). 

An n th order initial value problem for delay functional differ- 
ential equations is to find a function y in C b (t Q + a) such that 

y ^ n V(t ) = U(t, y ( * ) ) j t Q < t £ t Q + a. (7) 

Of course this type of problem is a special ease of the initial 
value problem (1); however in this form it is possible to give new 
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hypothesis on U so that (7) has at most one solution on any in- 
terval [t , t + e] for small £ > 0. The following theorem 
generalizes a result due to Wintner [8] which is given in ex- 
panded form in Hartman [4], page 34. 

THEOREM 3. Let h k (t , s), k = 0, 1, 2, * **, n-1 satisfy conditions 
(i) - (iii) of (4) and let U (t , x ( • ) ) satisfy 

n-1 f V a n s t . , 

|U(t, x ( • ) ) - U(t , y ( • ) ) | £ l X k (t) [ x (s) -y .( s ) |d n k (t , s ) 

k=o K Jt s K 

o 

for t in (t Q , t + a], where the n k (t) are non-negative functions 
such that 

n-1 ^Q +a (s-t yn-k 

J 0 x k (t) \ t cvkfr s) 1 1> t o <t 5 t o + a - 

o 

Then there is at most one solution to ( 7 ) on any interval 
[a Q , t Q + e] for small e > 0. 

The proof is an easy extension of the proof given in Hartman 
[4], page 558. 



10 


REFERENCES 


1. R. D. Driver. Existence and stability of solutions of a 
delay-differential system. Archives for Rational 
Mechanics and Analysis 10, 410-426 ("1962 ) . 


2. R. D. Driver. Existence theory for a delay-differential 
system. Contributions to Differential Equations 1, 
317-336 (1963). 


3, R. D. Driver. Existence and continuous dependence of 
solutions of a neutral functional— differential 
equation. Archives for Rational Mechanics and 
Analysis 197 149-166 (1965). 


4. Phillip Hartmen. Ordinary Differential Equations . 
New York, John Wiley and Sons, 19^4 . 


5. W. A. Coppel. Stability and Asymptotic Behavior of Differ- 
ential Equations. Boston, D. C. Heath and Company, 
1965. 


6. E. A. Coddington, and N. Levinson. Theory of Ordinary 

Differential Equations . New York, McGraw-Hill Book 
Company, 1955. 


7. F. Brauer, and S. Sternberg. Local uniqueness, existence 
in the large, and the convergence of Successive 
approximations. American Journal of Mathematics 80, 
421-430 (1958). 


8. A. Wintner, On the local uniqueness of the initial value 

problem of the differential equation d n x/dt n = f(t,x). 
Boll . Un. Mat. Ital. 11, 496-498 (1956). 



CHARACTERISTIC MULTIPLIERS FOR SOME 
PERIODIC DIFFERENTIAL EQUATIONS 

T. G . Proctor* 

1. Introduction. Let P(t) be a 2 x 2 matrix with elements 
which are continuous real valued functions of period T and 
consider the differential equation 

(1) x = P(t)x, 

where x is a vector with two components. It is well known [7] 
that there are numbers X 1 and X 2 , called characteristic multi- 
pliers, and corresponding solutions x-^(t), x 2 (t) of (1), called 
normal solutions, which satisfy for i = 1, 2, 

f T 

x^(t + T) = A. X .(t), -°° < t < 00 , X 1 X 2 = exp trace P(t)dt. 

If X 1 ^ X 2 and in some cases when X^ = X 2 any such normal 
solutions x 1 (t ) , x 2 (t) are independent. In this case a know- 
ledge of the characteristic multipliers and the values of 
x^(t'), x 2 (t) for 0 _< t _< T gives information for every solution 
of (1) for all t. It is clear that corresponding statements 
can be made for the second order equation 

~ \ + p(t) + q(t)y = 0, 
dt^ at 

*This work was supported by the National Aeronautics and 
Space Administration under Research Grant NGR 41-001-016. 



with continuous periodic coefficient functions p(t ), q (t ) 
since this results from the case 



y 



0 

1 

X = 

dt 

3 

P(t) = 

-q(t ) 

-P (t ) 


Calculation of the characteristic multipliers is not 
routine since in general one does not know even one non- 
trivial solution of (1), However it is possible to obtain 
convergent series representations for the solutions and thus 
calculate approximate values for the multipliers [2] , [8], 

An alternative procedure for obtaining the characteristic 
multipliers and the corresponding normal solutions for (1) is 
possible whenever an associated Riccati differential equation 

has a periodic solution. If we make the change of coordinates 
11 2 2 2 

x = z + y(t)z 3 x = z in (1) where y(t) is a solution of 
the Riccati equation 


||- = a(t) + b(t)y + c(t)y 2 , 


( 2 ) 


a(t ) = p 12 (t ) , b(t ) = p-, (t ) - p 00 (t), c(t ) = -p 01 (t ) , 


"11 


22 


21 


the differential equation in z can be integrated. This gives 

THEOREM 1. (a) If x(t) = column (x 1 (t), x 2 (t)) is a solution 

1 2 

of (1) then y(t) = x (t)/x (t) is a solution of (2) on any 

2 

interval on which x (t) does not vanish, (b) If y(t) is a 


solution of (2) on an interval I containing the number k then 
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x 1 (t ) 
(3) 

x 2 (t ) 


rt 

y(t) exp [p 21 (s)y(s) + p 22 (s ) ]ds 

■t 

exp Cp 21 (s)y(s) + p 2 2 ( s ) ] d s 

j k 


is a solution of (1) on I. (c) If y (t ) is a solution of 

(2) with period nT and f is the mean value of p 2 -^(t)y(t) + 

P 22 (t ) over the period nT, where n is a positive integer, 
riT f* 

then e is a characteristic multiplier for (1) for the 
period nT and (3) is a normal solution of (1) corresponding 
to this multiplier. 

In sections 2 and 3 of this paper we will prove two 
theorems with rather restrictive hypotheses which give the 
existence of a periodic solution of (2). The theorem in 
section 2 can be viewed as a special case of the theorem in 
section 3. We also mention other known techniques for 
constructing periodic solutions to the Riccati differential 
equation. 

Similar analysis for the differential equation (1) where 

P(t) is n x n and x is an n vector leads to the study of a 

matrix Riccati differential equation and the analysis is 

more difficult. The technique of using a Riccati differential 

equation has been used by Gelmand [ 3 ] and Andrianov [1] for the 

case of quasiperiodic coefficients p^(t) . 

2. b(t) has nonzero mean value. Let H be the set of all 

continuous real valued functions with period T, let b e H , 

ft 

let B(t) = b(s)ds and in this section we will assume 
J 0 

B (T ) ^ 0. Also suppose u, l e H and satisfy £(t) < u(t) for 
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all t and let K be the subset of H consisting of functions 
f where i,(t )<f (t )<u(t ), Further assume q(t, z) Is a con- 
tinuous real valued function defined for £(t) £ z £ u(t) 
such that for any fixed z, q e H and define CJ : K H by 


^h(t) = 


1 

-B (T ) 


-1 


t+T 


q ( s , h(s ) )exp 


b(v)dv ds. 


THEOREM 2. If (1) Ju(fc) £ u(t), Jl(t) £ 3l(t) for all t and 

d f m \ 

( 2 ) 


e x > 0 and z £ w implies q(t, z) £ q(t, w) or 


(3) e l < 0 and z < w Implies q(t, z) > q(t, w), then 


the sequences { J/u(t)} m=1 and { x(t)} m=1 converge 


uniformly to periodic solutions of 


(4) 


|f = b(t)z 


+ q(t. 


z) . 


Proof. The hypothesis (1) and (2) or (1) and (3) imply that 
if h 3 k e K and h(t) £ k(t) for all t then ^Th(t) £ J'k(t) for 
all t. Thus the sequences { (t ) }”_ Q and {^ m u(t)}”_ 0 are 
nondecreasing and nonincreasing respectively, uniformly 
bounded and equicontinuous . 

Remarks* If z(t) is a periodic solution of (4) where 
q (t j z) = c (t ) ( w ( t ) + z) 2 . 


w(t ) 


^BlfT 


-1 


rt+T 

Jt 


a(s) exp 


■t 

b(v)dv ds 
s 


we have that y(t) = w(t) + z(t) is a periodic solution of (2). 

If g(t) is a periodic solution of (4) where q(t, y) = 

2 

a(t) + c(t)y then y(t) is a periodic solution of (2). 
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Example . If we regard the right member of y = (y - 1) x 
(y - 2) as having period 2tt the corresponding transformation 
3 Is given by 


^h(t) 


-3t j*t +2 7T 2 

-fry , * 3s [h 2 (s) 

6 —1 t 


+ 2 ]ds . 


We note u(t) = 5/4 , S-(t') = 0 satisfy the conditions of the 
theorem and ) , ^f m u(t) converge to 1 as m 

Several theorems on periodic solutions of (2) are given 
In Hale [6, pp . 28-31] which result from the contraction 
mapping theorem and successive approximations. Using Theorem 
5.1 [6] a periodic solution exists if (1) |c(t)| is small 
enough (q(t , z) = c (t ) (w(t ) + z) 2 , y = w(t) + z (t ) ) or (2) if 
| a ( t ) | is small enough (q(t, y) = a(t) + c(t)y ). We state 
without proof an obvious modification (for the scalar case) 
of Theorem 5.2 [6] because of its relation to (1). Let 
h e H and have mean value zero. Let r and a be positive 
numbers and let q (t , y , y ) be a continuous real valued function 
for all t, -r £ y <_ r and |p| <_ c . Further let q be periodic 
for fixed (y, y ) and satisfy for |y|, |z| <_ p, |y| £ a 

|q(t, y, y) - q(t , z, y)| £ n(y, p)|y - z|, q(t, 0, 0) = 0, 

where p(y, p) is continuous and nondecreasing in y and in p 
for | y j < c, 0 £ p £ r , 



THEOREM 3* There is a 6, 0 < 6 < a, so that for |y| <6 


= y[b(t) + q(t, y , y) + h(t)] 


has a periodic solution which is the limit of successive 
approximations . 

COROLLARY. Let p(t, y), c(t, y) be continuous real valued 
functions defined on R x [0 3 a] which satisfy p(t, 0) = 
c(t, 0) = 0 3 p, c e H for fixed y, let q, r, h e H, h with 
zero mean value, r with nonzero mean value and let 


P(t) = P(t, 


S' q(t ) 
l-c(t , y) 


p(t, y) + h(t ) 


q ( t ) - r (t ) 


Then for y sufficiently small there is a periodic solution 
y(t) of (2) and e is a characteristic multiplier for (1) 
where f is the mean value of p 11 (t)y(t) + p 22 (t). 

By replacing t by t/y we see that an alternative state 
ment for this corollary is that there is a number X and a 
solution of 


x = 


q(t/y) 

-c (t/y, y) 


p (t / y , y) + h(t/y) 


q (t/y ) - r (t/y ) 


satisfying x(t + yT) = Ax(t) when y is small enough. 

Wasow [9] and Golomb [4] have developed a recursive 
scheme for constructing periodic solutions for the quasi- 
linear differential equation y = b(t)y + q(t, y, y) when 
b(t) has nonzero mean value and y is small. The resulting 
theorems provide another technique for obtaining periodic 
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solutions of (2) and thus obtaining characteristic multipliers 

for (1), Golomb [5] also shows that in certain cases this 

recursive scheme leads directly to a calculation of the 

characteristic multipliers and corresponding normal solutions. 

3* General Case. Let H / be any closed subspace of H in the 

1 f T 

uniform topology and for any f e H define Mf = f(t)dt. 

/ J ° 

Let £, u be in H with £(t) < u(t) for all t and let K be all 

functions in H y lying between £ and u. Further suppose 

q(t, y? u) is defined and continuous for £(t) < y < u(t) 

| v | £cr for some c > 0 and has period T for fixed (y, y). 

For any number a satisfying £(0) <_ a £ u(0) and | y | <_ a 

define 'I : K -»■ H by 

^ a, y '■ 


,f(t) = 


[q(s, f(s), y) - M{q(t , f.(t), y)}]ds. 


THEOREM 4. For fixed (a, y) if J (KjdH' 1 , J u(t) < u(t), 

3 . 3 . 3 y — 

£(t ) £ $ £(t) for all t and if y, z e K, y(t) < z(t) 

implies J y(t) < J z(t) for all t then the sequences 

d J 4 

} m _0 converge uniformly to functions 
£*(t, a, y), u*(t, a, y) satisfying 


u) = q(t, y(t, a, y), y) - M{q(t, y(t,a,y), y)}. 


The proof is identical to that of Theorem 2. 

COROLLARY. If £*(t, a, y) exists and M{q(t , £* (t ,a,y) , y)} = 0 
then £* is a periodic solution of 

ff- = q(t , y j y ) , 
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An analogous statement holds for u* (t , a, y ) . 

As an example of this theorem and corollary consider 
q (t , y, y) = y sin t, a = 1, and H / the set of all even 
functions in H, £(t ) = 1 and u(t) = e"*" -003 ^ . 

Hale [6, pp. 38—44 ] gives several theorems for periodic 
solutions of = b(t)y + q (t , y, y) in the ease where b (t ) 
has zero mean value and q is small for small y. These theorems 
again use successive approximation and the contraction 
mapping theorem. In the application of these theorems to 
the differential equations the symmetry properties of 
q(t, y, y) play an important role in the analysis to determine 
if the bifurcation equations M (t, y (t ) ) = 0 have a solution. 
One of the simplest cases arises when b(t) = 0 and q(t, y, y) 
is odd in t. In this particular case the recursive method 
developed by Wasow and Golomb also gives a convergent 
series expansion for the solution for small y. 

4. Final Remarks. Note that if T is the least period of 
P (t ) in (1) then the coefficient functions a, b 5 c in (2) 
have least period 0 or T/n for some positive integer n. 

By applying the construction procedures used in Theorems 
2, 3) or 4, we construct solutions of (2) with period T/n 
or 0. Then by replacing nT by T in Theorem 1, we obtain 
a characteristic multiplier for the period T . 
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For a given matrix Q(t) we notice that all matrices 
P(t) = Q(t) + K(t) have the same associated Rieatti 
differential equation (2) if K(t ) has the form 



If for some such K(t) equation (1) can be solved. Theorem 1 
may be used to solve equation (2). This in turn furnishes 
a method to solve x = Q(t)x. This remark applies even if 
none of the functions in Q(t) and K(t) are periodic. How- 
ever in the case of periodic P(t) suppose x = P(t)x has a 
normal solution, column (x (t), x (t)) where x (t) does not 
vanish. If X is the corresponding characteristic multiplier, 
y(t) = x (t)/x (t) and g = Mp(t), then Xe 6 is a characteristic 
multiplier for x = Q(t)x and a corresponding normal solution 
is given by (3) where p 22 is replaced by p 22 - p. 

Since the proof of Theorem la and b does not require 
the periodicity of P(t), theorems which imply the existence 
of bounded or almost periodic solutions of (2) or the 
asymptotic form of certain solutions of (2) also give 
information concerning the solutions of (1). 



10 


REFERENCES 


1. L. Ya. Adrianova, The reducibility of systems of n 

linear differential equations with quaslperiodie 
coefficients, Vestnik, Leningrad Univ , 17 ( 1962 ) , 
14-26 

2. N. P. Erugin, Linear Systems of Ordinary Differential 

Equations , Academic Press Inc., New York , 1966 . 

3. A. E. Gelmand, On the reducibility of a class of systems 

of differential equations with quaslperiodie coef- 
ficient^, Dokl . Akad. Nauk SSSR 116 (1957), 535-537. 

4. M. Golomb , Solution of certain nonautonomous dif fer- 

ential systems by series of exponential functions, 
Illinois J. Math. 2 (195977 45-65 . 

5. M. Golomb, Expansion and boundedness theorems for 

solutions of linear differential systems with 
periodic or almost periodic coefficients , Archive 
for Rational Mechanics and Analysis , 2 ("19 58), 
283 - 308 . 

6. J. K. Hale, Oscillations in Nonlinear Systems , McGraw- 

Hill Book Company, New York, 1963. 

7. P. Hartman, Ordinary Differential Equations , John Wiley 

and Sons, Inc., New York, 196 T! 

8. I. G, Malkin, Theory of stability of motion , AEC 

Translation 3352, U. S. Department of Commerce, 
Washington, 1952. 

9. W. Wasow, Solution of certain nonlinear differential 

equations by series of exponential functions, 
Illinois J. Math. 3 XI959), 45-65. 


CLEMSON UNIVERSITY 
Clemson, South Carolina 



T . G. Proctor , Periodic solutions for perturbed nonlinear 
differential equations . 

Abstract : The existence of periodic solutions of a 

periodically perturbed system of nonlinear differential 
equations is established. The construction of such solutions 
is proved in a more restricted situation. These results 
generalize well known results for perturbed linear differential 
equations. Examples are given. 



PERIODIC SOLUTIONS FOR PERTURBED 
NONLINEAR DIFFERENTIAL EQUATIONS 


T. G . Proctor* 

1. Introduction. In this paper we investigate the existence 

and construction of periodic solutions of a periodically perturbed 
system of nonlinear differential equations. The perturbed 
system is studied using an integral equation introduced by 
Alekseev [1], [2] which is a generalization of the variation 
of constants integral equation. The techniques used are 
analogous to those used in establishing the existence of 
periodic solutions in perturbed linear systems [33. Almost 
periodic perturbations of nonlinear systems have been studied 
by May [4] using a similar technique; however, our systems 
will not necessarily meet his requirements. 

2. Existence of periodic solutions. Let P, 6 , a be positive 

numbers with 6 < a, let and be the closed spheres of 

radii 6 and a respectively in R n and let f(t,x) be a C“ function 

n 

from RxS^ into R with period P in t. We make the following 
assumptions concerning the function f. 

(i) y in implies the solution x(t , t , y ) of the 
unperturbed initial value problem 

X = f(t,x), x(t) = Y 


Department of Mathematics, Clemson University, Clemson, 

S . C, 29631. This work was supported by the National 
Aeronautics and Space Administration under Research Grant 
NGR 41-001-016. 
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exists for 0 <_ t-x <_ P. 

(ii ) There is a set K C. so the function F given by 

F ( y ) = Y - x ( P , 0 , y ) 
is a homeomorphlsm of K onto F(K). 

(ill) f(t,0) = 0, 0 < t < P. 

We denote the continuous functions with period P from 
R into by Let the perturbation function g(t ,x , e ) be a 

continuous function from RxS 5 x[0,Eq] into R n and satisfy 

(iv) g has period P in t, g (t ,x, 0 ) = 0 . 

( v ) For k in S , and 
o 

P 

Y (k j e ) = |^-(P , s ,k(s) )g(s ,k(s ) , e )ds , 

J 6 Y 

0 

we assume y(k, e ) is in F(K) for 0 <_ e < 

Let S be the Banach space of all continuous functions with 
period P from R into R n with the supremum norm and let T be 
an operator from S g into S defined by 

t 

Tk(t) = X (t j 0 jF _ 1 ( y (k 9 e ) ) + |^-(t ,s ,k(s ) )g(s,k(s),e)ds 

J dY 

0 

for 0 <_ t <_ P, 0 • <_ c < Gq • 
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THEOREM 1. For e small enough T has a fixed point y (t ) . 

Further y(t) Is a periodic solution of the perturbed differential 
equation 

y = f (t ,y ) + g(t ,y , e) . 

Proof . The hypothesis on f implies x(t , t,y) and t , t , y) 
are continuous for y in S ^ and 0 <_ t-t <_ P . Conditions (iii ) 
and (iv) imply 


lim x ( t } 0 } F 1 ( y ( k , e ) ) ) = 0 
£->■0 

uniformly for 0 <_ t <_ P and k in so choose so that 
TS <c 5 for 0 < e < £ n . It is easily checked that T is 
continuous, is closed and convex and TS^ has a compact 
closure in S g ; therefore by the Schauder fixed point theorem 
there is a function y(t) = Ty (t ) . This representation shows 
y (t ) has a derivative which is given by 

y (t ) = f ( t , x ( t ) ) + g (t ,y (t ) , e) + 

t 

f 

f x (t , x (t , s ,y (s ) ) )-g£(t , s , y ( s ) )g(s ,y ( s ) , e)&s 

0 

where x(t ) = x(t , 0 ,F”^(Y(y , e) ) ) . Also we have 



f (t ,y(t )-f (t,x(t )) = 
t 


d_ 

ds 


0 

t 


•f (t ,x(t ,s ,y (s ) ) )ds 


.ax 


f x (t ,x(t,s ,y(s))^ r (t ,s ,y (s))[y(s)-f (s,y(s))]ds 


8y 


Therefore for 


W(t ) = y(t) - g(t,y(t))-f(t,y(t)). 


we have 


W(t ) = - 


0 


C\ -y 

f x (t ,x(t ,s ,y (s ) ) »s ,y (s) )W (s )ds 


The only solution of this equation on [0,P] is W(t) = 0. 

Example 1. The initial value problem x = x^, x(0) = y has 
solution 


x <ut,y> = 1- T (t-T) > -h - y - h 


±_ 0 < t-T < P, 


For K = [ -ip.,0] the function 


F ( Y ) = Y ~ 


1-yP 


Y in K 


is continuous and has continuous inverse 


a-/ a -4 a/P 


F 1 (a) = 
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for -1/6 P < a, <_ 0. The requirement (v) on g can be written as 

k in £ implies 
2P 

P 


-1/6P < 


— - — — ■ — 3 ' g (s ,k(s ) , e )ds 

^[l-k( s ) (P-s ) j 


< 0 . 


Theorem 1 establishes the existence of a periodic solution 
to a perturbed nonlinear differential equation; however, no 
construction is given for such a solution. In many eases it 
seems unrealistic to suppose that T is a contraction operator 
since may not be Lipshitz as in the example. The following 

section does provide a method to construct periodic solutions 
in a special situation. 


3. Construction of periodic solutions. Let ft be a region 
in R n , and let f(t,x) be a C ' function from Rxfi into R n with 
period P > 0 in t and let £(t), u(t) be continuous functions 
from R into ft with period T with £(t) <_ u(t), 0 <_ t _< P where 
a vector inequality l <_ u means the components u^ of the 

vectors l, u satisfy l ^ _< u^, i = 1,2, .. . ,n. Let 

S = {x in R n : £(t)<x<u(t) for some t) 


then for y in S we assume 
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( i ' ) the solution x(t ,x ,y ) of the initial value problem 

x = f(t,x), x ( x ) = y 

exists for 0 <_ t-x f_P, t any real number; and 

( ii ' ) the function F given by 

F(y) = y - x(P,0,y) 

is a homeomorphism of S onto F(S). 

Let g(t,x) be a continuous function from RxS into R n with 
period P in t and let S* be the set of continuous functions 
k(t) from R into R n with period P satisfying £( t) <_ k(t) <_ u(t) 
for all t . For k in S* and 

P 

y (k ) = |y(P,s,k(s))g(s,k(s))ds 

0 

we assume y(k) is in F(S). And we define an operator T on 
S* into 5 by 

t 

Tk(t) = x(t ,0,F -1 (y(k) ) + |^-(t , s,k (s) )g(s ,k( s) )ds 

0 


for 0 < t < P. 
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THEOREM 2. If l(t ) <_ T£(t), Tu(t) <_ u(t) for 0 < t <_ P and if 
k, h in S* with k(t) <_ h(t), 0 <_ t £ P implies Tk(t) <_ Th (t ) , 

0 £ t <_ P, the sequences {T m £} m ~Q, T m u converge uniformly 

to fixed points of T. If y is a fixed point of T then y(t) 
satisfies 


y = f (t ,y) + g(t ,y ) . 

Proof. The sequences (T m £(t)} °° n and {T m u(t)} ” are 

m= u m= U 

nondecreasing and nonincreasing respectively, uniformly bounded 
and equicontinuous . Hence they converge uniformly to limit 
functions which are fixed points of T. The proof that such a fixed 
point is a solution of the perturbed differential equation is 
identical to that given for Theorem 1. 

Example 2. The initial value problem x = x(l-x), x(0) = y has 
solution 


x(t,x,y) = — ye - 

1+ Y (e t -1) 


for y > 0, 0 < t-x < P. For y* = (e^^-1 )/(e^-l) 


the function 


F(y) = y H5 , y* < y £ 1, 

l+y(e -1) 
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and let £(t) = 0, u(t) = 1 , 0 <_ t <_ P. It is easy to verify 
the conditions of the theorem if P is small enough. 


4. Final remarks. The hypothesis ii and ii ' of Theorems 1 
and 3 respectively is analogous to the noncriticality requirement 
[32 made for unperturbed linear systems. The existence of 
periodic solutions to perturbed nonlinear systems corresponding 
to the critical case can be treated using the methods above. 
However, if one imposes hypotheses similar to those used in 
Theorem 1 the question of a solution to the resulting bifurcation 
equations is not easily established. This follows since we 
do not know the dependence of the fixed point y(t) as a function 
of e. Hypotheses similar to those used in Theorem 2 are extremely 
hard to verify since this requires that the difference between a 
function g (t , x ) and its mean value be monotone over a class of 
periodic functions x(t). 
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A FAMILY OF SOLUTIONS OF CERTAIN 
NONAUTONOMOUS DIFFERENTIAL EQUATIONS BY SERIES 
OF EXPONENTIAL FUNCTIONS 

by T. G. Proctor and H, H, Suber 
1. INTRODUCTION 

We consider in this paper the construction of solutions 
for certain nonautonomous differential equations. The first 
result makes use of a technique developed by Golomb [ 5 ] 
and Wasow [10] for constructing solutions of some non- 
linear differential equations by means of series of 
exponential functions. The technique as employed here 
gives explicit formulae for' a family of periodic solutions 
of a Ricatti equation with odd periodic coefficient and 
finite Fourier series expansion. 

Following this is a theorem concerning the existence 
of a family of almost periodic solutions of the vector 
differential equation 

y = g(t, y). 

Here y is an m-vector; g(t, y) is quasi-periodic and odd in 
t and satisfies certain other conditions. (A quasi-periodic 
function is a function almost periodic in t with a finite 
base of frequencies w-j Wg , • * * , u> n . ) The theorem is 
a generalization of a result concerning periodic solutions 
when g(t, y) is periodic in t [2], in particular, the 
Ricatti case mentioned above. The proof of the theorem 
utilizes a method devised by Kolmogorov [7] to overcome the 



problem of arbitrarily small divisors and gives a method of 
constructing approximations to the almost periodic solutions. 
Since we assume that g has a finite base of frequencies we 
can present the system of equations in an autonomous form by 
considering a higher dimensional version of the differential 
equation. The theorem is as follows: 

Let x be an n vector, let y be an m vector and consider 
the differential equations 

x = w , y = f(x, y) , (1.1) 

where f(x, y) = -f(-x, y), and where the components of 
f(x, y) are analytic for |y| £ R^ and j Imx | <_ R 2 , and 
where f has period 2ir in each of the components of x. 

Suppose that the vector ft = (u^, u> 2 , • • •, w n ) satisfies 
an inequality 

!k-n| > —5- (1.2) 

- |k| v 

for some positive constants K and v and all vectors 
k = (k^, k 2 , k n ) with integer components where 

n n 

k*ft = l k.m. and Jk| = T |k.| ^ 0. Then if f(x, y) 

1=1 11 i=l 1 

is sufficiently small for |y| £ and | Imx | £ R 2 , there 
is a neighborhood of y = 0 such that all solutions of (1.1) 
starting in this neighborhood are almost periodic with base 
frequencies ft. 
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We note the requirement on ^ (inequality (1.2)) is not 
stringent. If v > n, such a constant K exists for almost 
all ^ (in the sense of Lebesgue measure) [2]. The proof 
of this theorem is given by constructing an infinite se- 
quence of coordinate transformations so that in the limiting 
set of coordinates the differential equations can be in- 
tegrated in a neighborhood of y = 0. 

The last section gives an application of the results 
mentioned above. 


2. THE PERIODIC CASE 

Let 1^ denote the set of all k-tuples of non-negative 

integers. The elements of 1^ can be counted according to 

the following technique. Let (n-,, n 2 , • • *, n, ) e I, , let 
n l n 2 n i n k 

m = 2 3 * * • Pj_ • • * p^ where p^ is the ith prime 

number i = 1, 2, • • • , k, then denote (n^, n 2 , • • •, n^) 
by N m . Note that N a = (0, 0, • • •, 0), N 2 = (1, 0, • • • , 0), 
etc. The natural ordering of the non-negative integers then 
orders all the elements of 1., N^, N 2 , * * v. For 
N n s ^m e I k we make the following remarks: 


1) 

m is prime iff N i 

m 

s of the form (0,0, 

• * • ,i, ••• 0) 

il) 

Define addition in 

1^ component-wise. 

Then 


N + N = N 0 iff nm 

= A. 


iii ) 

We say that N n < N m 

iff n < m and N < 
— n 

N 

m 


iff n < m. From ii ) it is clear that 

N + N = N „ implies that N < N„ and N < N. . 
n m Z n — Z m — Z 


3 



For a) real let P, = 

(oj, - to, 2o)j 

3 

C\J 

1 

• 

, km, -km). Now 

for N m = 
m 

: (n 1# n 2 , * 

• • “n 

5 2k-l’ 

n 2k ) 

£ X 2k 

let m'' be the 

integer 

so that N ^ 

“ ( n 2 3 5 

• • 

5 n 2k’ 

n 2k-l ) * We 

observe 

that N = 

m 

-N Let 

m 

°k 

denote 

the class of 


odd periodic functions of the real variable t with period 
2irw which have finite Fourier series containing only terms 

of the type f.e J for j = ±1, ±2, • • • , ±k. Using 

3 

the above notation we may represent functions in 0^ in the 
form 

iN • fit 

f(t) = l f e (2.1) 

N el„, 
n 2k 

n prime 

where f = -f since f is odd. 
n n 

Theorem 2.1 . For n > 0 sufficiently small the differ- 
ential equation 

y " = n(a(t) + b(t)y + c(t)y 2 ), (2.2) 

where a, b and c e 0^, has an even periodic solution y(t) 
of the form 

iN «S2t 

n • (2.3) 

Proof . Assume that (2.2) has a solution of the form 
indicated. Then formally we have 


y(t) = 


l 

N n £l 2k 


y e 
J n 



N el 0 , 
n 2k 


iN *fiy e 
n J n 


IN *Qt 
n 


n l 


N n £l 2k 


a + 7 yb+ y y n y c 

n L J i p „ A J £ J m p 

mp=n £mp=n * 


p prime 


p prime 


iN • fit 
n 


(2.4) 


where a., b . and c . , j prime are the coefficients of the 
J J J 

series representations (2.1) of a, b, and e respectively, 

Now in case N *fi ^ 0 we write 
n 


y 




a. 


iN *fi a n 5 


n prime 


n 




(2.5) 


n 


iN *fi 
n 


I y m b D + I yoy m c D 

mp=n m p £mp=n * m p 
p prime p prime 


n not prime 


and for N n *fi = 0, 


y n = 0 . (2.5a) 

iN *fit 

Suppose that the terms containing e n on the right 
side of (2.4) vanish whenever N n *fi = 0. Then by remark iii) 
above we see that (2.5) defines y n recursively so that (2.3) 
will be a formal solution. To show that this is indeed the 
case we present the following. 

Lemma 2.2 . Let {y n } he the sequence of numbers defined 
by 2.5. Then y n = y n ^. 
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Proof . The proof is by induction. If N n *fi = 0 then 
clearly y n = y ^ and so in particular for n = 1 . Now 
suppose that y m = y ^ f° r all m < n. Then we may write 


n 


iN -a 
n 


T y b 

L J m. p 
mp=n ^ 


I 

£mp=n 


y £ y_ c 


m p 


-iN n ^ *Q 


I y m^ ( ~ b n^ + l y ^ y m^- C n^ 


[mp=n 


P' ' » " A m' ' "p- 

Amp=n 


where p is prime. But mp = n iff iti'p" = n' . So we see 
that y n = y n > which proves the lemma. 

Lemma 2,3 - If N n e I^ k is such that N n *fi = 0 then 

I y m b p + l y £ y m C D + ^ . y m b n + l . y A y m C p = °» 

mp=n m p Amp=n * m p mp=n m p Amp=n" * m p 


where p is prime. 


Proof . By Lemma 2.2 we have y n = Y n ' for all n and 
since b and c e 0^, we may write 


l y m b D = I = - l y m b D 

mp=n m p mp=n m p mp=n" m p 


and 


trnp- 


I y £ y m c n = I y Jl" y m" (-C D" ) = “ l y £ y m°n’ 
=n * m p Amp=n * m p Amp=n" % m p 


Amp = 


where p is prime. 
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Now in order to show that (2,2) is a solution of 
equation (2.1) we will prove that formal series (2.3) with 
y n defined by (2.5) converges uniformly and absolutely for 
all t and n sufficiently small. 

Let C represent the complex plane, for 
z = (z^ , z 2 , • • • , z 2k ) e C 21c and N n e I 2k define 

Z j 

2 k d 

mapping G into C is analytic in the poly disk 

{z e C 2k : |z i | £ r } of radius r about the origin iff f 

has the representation. 


A function f 


_ n _ “1 “2 

z — z z 2 


2k 

z 2k and z = max 


f(z) = I 


N 


N n eI 2k 


a n z 


n 


where the sum is uniformly and absolutely convergent in 
the polydisk. In case f is analytic, Cauchy’s inequality 
gives for |z| <6 


N 


a n | < M/6 


n 


where |N n | = n^ + n 2 + 


Now for z e C 2k let 


+ n 2k and M 


sup | f ( z )[, 
|z|=6 
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Z . 

J 


3 


a* ( z ) 


2k 


3=1 


|a j 


2k 

b*(z ) = l | b . | z , 
3=1 J J 


2k 


(z) = l | c . | z . ; 


3=1 


J J 


and let u(z ) = f (a* (z ) , b*(z), c* (z ) ) where 


f (a,b,c) = n 


1-rib l-r)b 


2ric 2nc 

na 

l-nb 5 


1 - 


4n ,2 ae 


- :J s 


(l-nb )‘ 


3 c o. 


c = 0 . 


( 2 . 6 ) 


Note that u Is the solution of the equation 

u(z ) = n[c*u 2 (z ) + b*u(z) + a*] (2.7) 

which vanishes when a* = b* = c* = 0. We see that u(z) 

Is an analytic function of z in any region which does not 
include zeros of the function 

g(z) = (l-nb* (z ) ) 2 - 4n 2 a* (z )c* (z ) . 

Now for 6 > 0 choose n Q > 0 so that | g(z ) | > 0 
whenever |z| < 1 + 6; e.g. for L = max{ | a . | , | b . | , | c . | } , 

3 3 3 3 

let n o < BLk(l+6) » then %l a * I > 1 5 n ol c *l < W and 
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we see that in this case |g(z)| > 0. Now for all n, 

0 < h < n Q we have u analytic in the polydisk |z| <_ 1 + 6. 
Hence, in this polydisk u has the representation 


u(z) = l u z n , 
N n eI 2k 


where 


I u n | < M/(l+6) 


with M = sup | u ( z ) I . 

| z. | =1+6 

On the other hand, substituting frtim (2.8) into (2.7) 
and using (2.6) we obtain 


nl a n l 


n prime 


n I -u m |b I + n l u U c 

mp-n m p £mp=n * m p 


n not prime , 


p prime 


p prime 


Comparing this with the recursion formula (2.4), with 
y-^ = 0, we see immediately that 


i y nl 1 ml u nl > n = 1 r 2 > ' * * ‘ 

IN • fit 

Since I e 1=1 for all t we have 


9 



3 


N 


l 

n e ^2k 


y e 


IN -fit 
n 


< 1 
— w 


U 


N n £l 2k 


n 


< M T 1 

w N el TC~ 

W n X 2k (1+6) n 


^ M (l+6'i 
- w *> 6 ' 


2 k 


which not only proves absolute and uniform convergence, 
but also gives a bound for the solution y (t ) . 

Remarks : 

i) In the proof of Lemma 2,2 we showed that 

y n = Y n ' for all N n e such that 

N ¥■ 0. Prom this we conclude that 
n 

the solution found above is even in t. 

Note also that the solution has zero 
mean value . 

ii) The particular order relation used here for 
1^ is not essential to the proof. See 
Golomb [5] and Wasow [10 ] for different 
schemes . 

iii) It is possible to use the result in this section 

directly to obtain solutions with mean value other 
than zero. Let f(t, y) represent the right side 
of equation (2.1) and suppose that y(t) is the 
solution of 
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y' = .f(t, y) 

given above.. For any fixed constant c, let 
z = y + c in (2.1). The theorem gives a 
technique for obtaining a solution of the 
new equation 

2" = f*(t, Z) , 

where f*(t, z) = f(t 3 z - c) with zero mean 
value. This in turn gives a solution to the 
original equation with mean value -c . 

iv) Let y be an n-vector, let p (t) be an n-vector 
with components p^^'^t) e 0^, j «= l 3 2 3 *** 3 n 3 
Z - 1 3 2 3 ,, ‘. Then the differential equation 

N. 

y' - n l p„(t)y Z (2.10) 

» . . 

where the right side converges for |y| < r may 
be solved using the techniques of this section. 

The only essential difference occurs when one 
attempts to find an analytic solution of the 
corresponding equation (2.8). Here one may use 
the implicit function theorem to show existence 
of such a solution for n sufficiently small* 

v) The existence of periodic solutions of equations of 
the form (2.10), for n small is shown by Hale 
[6, p. 45]. 
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3. THE QUASI PERIODIC CASE 

For any positive integer n let denote the set of 

all n-tuples of integers, for a = (a^, ctg, ' * *, a ) e J 
n 

let | a j = £ | a. j and let C m be all m vectors (y 1 ,*** J y ) 

i=l 1 i m 

where each component is a complex number. For simplicity we 
will treat only the case where y is a m = 1 vector. 

We shall be concerned in this section with functions 

n 1 

defined and analytic on (x, y) subsets of C x C into 
C 1 which are periodic of period 2 tt in each component of x. 
These subsets will be of the form 

D(r, p) = { ( x , y ) e C 1+n : j Imx | £ p, jy| <_ r] 


where the norm j | of the vector x denotes the maximum 
of the absolute value of its components. We denote the class 
of such functions by P(r , p) and note that any 

g e P(r, p) has a Fourier-Taylor series representation 


gU, y) 


M?|$|=o feaf3 


ia*x 3 
e y 


3 


where the g a ^ are complex numbers and where the sum is taken 
over all a e J n and 0 e 1-^. 

Several lemmas are listed below without proof. The 
proofs are elementary and are similar to those given in [2 ] . 
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Lemma 3 • 1 • Let h e P(FL , H ? ) and let jh(x,y ) | £ M, 
M > 0 In D(R^, Rg ) • The Fourier -Taylor coefficients , given 
by 


h 


ag 


1 

( 2 TT ) n 


hg(x)e ia x dx-^dx 2 • ’ *dx n , a e J 


n- 


3 e I 


where 


h a ( x ) = 


1 

3! 


9y 


3 b (x, y) 


y=o 


and where the jth Integral is taken from x^. 
satisfy the inequality 


0 to x. 
J 


2tt, 


I h 


Me 


- | a |R2 


a 3 


Rn 


If h(-x 




have h 



, y ) = -h(x, y) in the above we have 

and conversely. If h(-x, y) = h(x, y) we 

h Q and conversely, 

ag 


Lemma 3 » 2 . If the elements of the sequence 
{ h a g } a e J n , 3 e I 1 satisfy 


- a R, 


h 0 \ < 

1 a 3 1 — 


Me 


R-, 


9 


then 



h(x,y ) = 


a , 3 =0 


h „e 
a3 


iavx 3 

y 


— R 

is analytic for |y| <_ R^e , | Imx | <_ Rp - <$ for any positive 
6 < 1 such that Rp - 6 > 0; and in this domain we have 


h(x. 


y) 1 



2 » +1 m 

s n+1 


Lemma 3 . 3 ■ For all positive numbers m, v, 6 we have 


m 


v 



e m« 

6 V 


Lemma 3.4 . (Cauchy’s Inequality). If the complex 
valued function h(z ) is analytic and bounded by M for 
| z | £ R, M ,R > 0, then for |z| <_ Re” , 0 < 6 < 1, we have 

dh 2M 
dz - R5 

i 

The proof of the main result of this section depends 
almost entirely on the following considerations. 

Let f e P( Rj > R 2 ^ and satisfy f(-x, y) = -f(x, y ) , 
let f satisfy 

| f ( x j y) | < M = 6 2(n+v)+7 (3.1) 

in D(R^, Rp) where 6 is specified below, let w satisfy (1.2) 
for some positive constants K and v and all a e J n with 
| a | ^ 0 and consider the differential equations 

x = a), y = f (x, y) . (3.2) 


Ik 



2 5 there exists 


Lemma 3 • 5 • For each x in | Imx | <_ R 2 - 
an invertible transformation U defined on a subset of C 
into C, given by U ( r> ) = y where 


y 


n + u(x, n), |u(x, n ) j < ~ n+v+2 

5 

.-46 


(3.3) 


for | n | £ R-^e -40 , | Imx | <_ Rg - 26. Putting ri = U“^(y) 
in (3*3) we obtain the differential equations (3.2) in the 
new coordinates 


x = w , 


n * f *(x, n); 


( 3 -^) 


-4 6 


and in D (R^e , R 2 - 26) 

|f*(x, n) I < m 3/2 . 

Further we have f*(x, n) = -f*(-x, q) and 


(3.5) 


ii (S 

f* e P(R^e~ , R 2 - 26), Here 6 is taken as a positive 


number satisfying 

l/n+v+4 

6 < mini 



R., K „ 

r 1 r® i v K 
3 n 2 3 2n+3 l v J 3 


2n+l '-v 


_ v l/n+v+6 

(!) . (?) 


Proof , a) Definition of u(x, n ) : Choose u(x, n ) as 

the solution with mean value zero of 


9u 

3x 


w = f ( x } n) 


(3.6) 


where is the vector with elements in the ith row and 


jth column. This gives 
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U(x, r) = I 

I a 1 , ] 6 j =0 “ B 


where 


u 


ag 

a$ ia*w 


a * °> u o B " 0 


Since |f(x, n)| £M for |y| £ R -^, | Imx | £ R 2 , Lemma 3.1, 
inequality (1.2) and the above imply 

-Mr, 


u „ < 

rv R 


Me 


a 


v 


a$ 1 - R ]0 


K 


By Lemma 3 • 3 


- | a j (R 2 -6) 

i i ^ M r v^ v e 

I U mR I — (eJ 


ae ' " K6 V ^ 6 ‘ 




3 


Hence, using Lemma 3.2 we have u(x, n) defined and analytic 
6 

for | n | £ R^e - , | Imx | £ R 2 - 2 6 and bounded in this domain by 


u (x, n) | £ 


,2n+l 


t n+l+v K 


M (|) 


v 


(3.7) 


Thus, inequality (3.3) is valid and we note that 
u(x, n) = u(-x, n) and u e P(R^e , R 2 - 26). 

b) The transformation U : By equation (3.3) the 

p j> 

set D = {q e C : |n| £ R^e - } is mapped into a set 
containing A = {y e C : ]y| £R 1 e - ^ < ^}; i.e. U(D) D A. 
Since 


l6 



for |ri | ^ R^e j I Xmx | £ Rg - 26 , we see that U 1 is 
defined on A. Thus for |n| £ R-^e , | Imx | £ R£ - 26, 
u(x, n) is defined and (3.7) holds. 

c) The function f#(x, n): Substituting from (3-3) 

into (3*2) gives 

(1 + |^-)n = f (x, y) - f(x, n), 

so that 

f*(x,n) = (l + |^ (x,n) ) _1 (f (x,n + u(x,n) ) - f(x,n))j (3.8) 

and we note here that -f*(-x,n) ” f* (x,n) • 

Now we have 


| (1 + 


3u 

3n 


)| 


-1 


< 2 


(3.9) 


hence 


|f*(x, n) | £ 2 |f (x,n + u ( x , n) ) - f(x, n) I • 

But 

r* Q 2n+1 .Ur v 

| f(x,n + u(x,n) ) - f (x,n)| £ sup { lff~l I (5) * 
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<5 

where the supremum is taken over jy| £ R^e , | Imx | £ Rg • 

By Cauchy’s inequality 


3f 2M 

3y - R 1 6 ’ 

thus 


f*(x. 


n) < 


KR 1 6 


2n+3 

n+v+2 





and the proof of the lemma is complete. 


Theorem 3.1 . Let f be as in Lemma 3-5. Then if M 

(and thus 6) is sufficiently small for each x in j Imx | £ Rp/2 

there exists an invertible transformation, V, defined and 

-R 2 

analytic on {n e C ; |nj £ R^e } into C m , given by 
V(n) = y where 

y = n + v ( x , n) • (3.10) 

Denoting the inverse transformation n = V - '*'(y) we obtain 
the differential equations (3.3) in new coordinates 


x = w, q = 0 . 

Furthermore we have v(-x, ri) = v(x, q). 

Proof. Choose 6, > 0 so that for 6. 

■ 0 



we have 
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I s • 1 


1=1 J 


R 2 

TT 


6 ll 



l/n+v+4 


R r K 


(!) 


K 


’ 2 2n+ 3 l v 


5 2 2n+l W 


(I) > ( 4 ) 


l/n+v+5 


Apply Lemma 3.5 Iteratively j times. Let (x, y) denote 

the function in the transformation of coordinates at the 

ith step and let f^(x, n ) denote the corresponding right 

side of the differential equation. We obtain the composite 

map F.(x 5 n) = y where 
J 

f. (x,n) = n + u. (x,n) + u._ 1 (x,u. (x,h) ) + • • •+ u-i ( x j n + u. (x,n) + * • • + 

u 2 ( x , n + u. (x, rj) + + u ^ ( x j n) ) , 

defined for |n| £ R-, exp [-4 > 6.], jlmxj < R„ - 2 ^ 6., 

1 1=1 1 ' 1=1 1 

where in the associated differential equations 


x - a) , n = fj ( x , n ) , 


the functions f . (x, n) satisfy 

J 


f, (x, n) I < m. = {. 2(n+v)+7 

J J J 


We observe that the composite transformations are defined 

-Rp 

for all j in | rj | £ P^e , | Imx | £ Rg/2 , and that 


I l u il < 

1=1 1 

thus the limiting composite transformation 
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(3.11) 


F ( x , n) = 11m F. (x, n) = n + v(x, n), 

3 + 00 J 

will exist in the above domain. In the coordinates defined 
by (3.11) the differential equation (3.3) becomes 

x = a) , n = 0 . 

4 . APPLICATION 

Adrianov [1] and Gelmand [4] outlined a procedure for 
finding a transformation x = Z(t)y so a given differential 
equation 

§f = Q(t)x, 

where Z(t) Q(t ) are almost periodic n x n matrices and where 
P satisfies certain conditions and x in an n-vector, becomes 

y = Ay, A constant 

in the new coordinates. We shall follow this procedure and 
use Theorem 3-6 to effect the same transformation in cir- 
cumstances where the earlier work fails to apply. 

Let H be the class of all functions 

f(t) = g(oo 1 t, ft) 2 1, ' * *, m n t) 

where g(u^, u 2 , * * * u ) is real analytic and has 
period 2 tt in each u^, i = 1, 2, • • • , n. Consider the 
differential equations 
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(4.1) 


dx^ 

= [q(t) + nq n (t)]x 1 + nq 12 (t)x 2 . 


dx 2 

= nq 21 (t)x ] _ + [q (t ) + nq 22 (t )]x 2 , 


where n > 0 and q, q. . e H and q. . 

J 1 J 

and where w = (w^, w 2 , • • • , w ) 

We make the change of coordinates 


(t ) = -q ±j . (-t), i j j = 1, 2, 
satisfies inequality (1.2). 


Xn = y n + Ty„, x 0 = y. 


(4.2) 


where t is any almost periodic solution of the differential 
equation 


T = n[q 12 + (q lx - q 22 )r - q 21 x ] 


(4.3) 


In the new coordinates (4.1) becomes 


dy l 

dF~ 


dy 2 

dt 


[q(t ) + n(q 11 - 

nq 21 y l + n(q 21 T 


T q 21 ) ]y x 

+ q 2 2 q)y 2 


(4.4) 


Theorem 3.6 guarantees that for n small enough equation (4.3) 
has almost periodic solutions, which belong to H. Equation 
(4,4) may now be integrated to obtain 

y 1 = a*(t )e ° c ± 

at b t 

y 2 = b* (t )e c 1 + b**(t)e c 2 
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where a Q and b Q are the mean values of q + nq-Q - nxq 21 
and q + nq 2 2 + hq 2 ]_ T respectively and a*, b* and b** e H. 

Reversing the change of coordinates (4.2) we obtain a 
fundamental matrix solution of (4.1) of the form 


<|>(t) = P(t)e At , 



5 


where the elements of P(t) are almost periodic and belong 
to H. 

The change of coordinates 


x = P(t )z 

in (4.1) yields 

Z = AZ . 

The case treated by Gelmand [4] required that the linear 
term in the resulting differential equation for x have mean 
value which dominates the other elements in order that there 
exists an almost periodic solution of this equation. Thus 
Theorem 3.6 permits us to consider a new situation. 

If in 4.1 we require that q, q^. e 0^ (defined in 
section 2) i, j = 1, 2, then we may use Theorem 2.1 to obtain 
an explicit representation for periodic solutions of equation 
(4.3) for r\ sufficiently small. Then we may integrate equa- 
tions (4.4) and obtain explicit solutions of (4.1). Note 
that if 
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T (t ) 


IN 


n= 1 


T n e 


n 


fit 


is the solution of (4.3) given by Theorem 2.1 and t 1 = 0 
we have periodic solution of (4.1). The existence of these 
solutions was shown by Epstein [3] • 
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INTRODUCTION 


This dissertation establishes the existence of integral 
manifolds for a system of perturbed nonlinear differential 
equations in a neighborhood of a critical point. 

The notion of an integral manifold for a system of 
differential equations, as used in this discussion, is 
given here. Let U be an open set in R p and consider the 
system of differential equations 

(1) z = Z ( t , z ) , 

where z, Z are real p-vectors and Z is continuous in RxU . 

Let z(t,r ,c) denote solutions of (1) satisfying the 
initial condition z(x,x,t) = c. For a fixed number t 
suppose there exists a set S + which can be represented in 
the form 

S + = { ( t , z ) in RxR p | z 1 =v i (t ,z k+1 , . . . ,z p ) , i=l,...,k, t^x}, 

where v = (v-^,...,v k ), k < p, is defined and continuous 
on [x,»)xV, V a subset of The set S + will be 

called a positive (integral) manifold for system (1) if 
any solution, z(t,x,s) of (1), with (x,c) in S + exists 
for all t >_ x and is such that (t , z (t , x , s ) ) is in S + for 


all t > x. 
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Similarly a set 

S" = { ( t , z ) in RxR p | z£=tu(t 3 z k+ ]_> • » • > z p) > i=l , t<_t}, 
where u = (u,,...,u ), l < p, is defined and continuous 

■ X/ 

on (-°°,x]xW, W a subset of R y , will be called a negative 
( integral ) manifold for system (1) if any solution, 
z(t,T,c) of (1), with (x,£) in S~ exists for all t £ x 
and is such that (t ,z (t , x , ? ) ) Is in S _ for all t £ x. 

Suppose that system (1) possesses a positive manifold 
S + with v bounded on its domain. Consider the perturbed 
system 

(2) z = Z(t,z) + Z*(t,z), 

where Z* is a real p-vector. Sufficient conditions on 
the function Z* are given to insure the existence of a 
positive manifold 

S* = { ( t , z ) in R 1+p | z ± =v* (t ,z k+1 , . . . ,z p ) , i=l, . . . ,k,t>x}, 

with v* = (v*,. . > ,v*) bounded on its domain. A similar 
result is obtained for negative manifolds. 

A restatement of the problem in a special case 
serves as a model. The system 
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3 * 3 

x * -x-\ y = y , 


has local solutions 


x ( t , t , a ) = 


9 1/9’ y(t } T,0) o 1/9 

(l+2a (fr-T)) 17 ■ (1-23^ (t-x) ) 17 ^ 


The phase portrait manifold of (3) is given in Figure la 
Notice that 


S + = { (t ,x,y ) | 0<x, y=0 , t >0} , 

is a positive manifold for (3) and 

S - = { (t ,x,y ) | x=0 , 0 <y , t<0}, 

is a negative manifold for (3). The problem here is to 
find sufficient conditions on the perturbation functions 
X(x,y), Y(x,y ) so that the system 

x = -x 3 + X(x,y), y = y 3 + Y(x,y), 


has positive and negative manifolds 
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S # = { (t , x , y ) | 0<_x<e , y - v + ( x ) , t > 0} , 

S; = { (t,x,y) | x=v~(y ) , 0£y<n , t<_0} , 

with v + and v” bounded on their respective domains. The 
phase portrait of this system may then be represented 
schematically by Figure lb. 



In previous work on this problem the typical approach 
has been to assume that system (1) has a manifold S. In 
many cases solutions on S are periodic in t . Then under 
the proper conditions, see [9], there exists a reversible 
transformation, T, defined in a neighborhood of S which 
carries (1) into a system of the form 



y = By + G (t , 0 ,x,y ) , 


6 



where d is a constant ^-vector , 0 an Ji-vector , x an 
m-vector , y an n-vect or , £+m+n = p. Here A is an mxm 
matrix with eigenvalues which have negative real parts, 

B is an nxn matrix, with eigenvalues which have positive 
real parts, and the functions 0, F and G satisfy certain 
smoothness and order conditions and are multiply periodic 
in 0 with period w = on > o, 1 = 1,2,.... 

i.e. 0(t,0+w,x,y) =* 0 (t , 0 ,x,y ) where 

0 + U3 = (0^+10^, . . . ,0^ + W^) . 

Equations (4) represent system (1) in a neighborhood 
of a critical point, periodic orbit or periodic surface 
depending on whether the 0-equation is absent, i = 1, 
or H > 1, In case 1 = 0 system (4) may be considered a 
linearization of the original system. The transformation 
T sends the perturbation function Z* into the triple 
(0# ,F* ,G* ) and the problem then becomes one of finding a 
manifold with some specified properties for the system 

0 = d + 0 ( t , 0 , x , y ) , 
x = Ax + F (t , 0,x,y) , 
y = By + G(t , 0,x,y) , 


( 5 ) 



7 


where 0=0+0*, etc. One usually begins with a system 
in the form of (5). Functions describing the manifolds 
for (5) are then obtained as solutions of a certain 
improper integral equation, see for example ([13], p.137), 
formulated by the use of the classical variation of 
constants technique. The linear term in (5) gives rise 
to factors of the form exp(tA) which guarantee the 
convergence of these integrals. (An exception to this 
is work by Kelley [153 ) • This approach was introduced 
by Bogolinbov and Mitropolsky [3] s [4] and subsequently 
developed by Hale [12] and Kelley[15]. Of interest here 
is the case in which the linear terms in (5) are replaced 
by nonlinear terms. Again one is led to consider an 
improper integral equation, obtained in this case from 
a generalization of the variation of constants formula 
due to V. M. Alekseev [1], Special hypotheses are 
Introduced to insure convergence of the improper 
integrals . 

Some work on the nonlinear problem has been done by 
L. E. May [173. The hypotheses used by May seem to admit 
only systems which behave in an essentially linear manner 
in the unperturbed state. The theorem presented here, 
which includes the model problem above, requires no 
linear or approximate linear behavior in the unperturbed 


state . 
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Other people who have worked on the problem of 
integral manifolds are Diliberto [10] , [11], Levinson 
[ 16 ] and Sacker [ 1 8 ] . 



PRELIMINARY CONSIDERATIONS 


The symbol || will be used to represent any fixed 
norm on R^ , R m or R n as well as the corresponding 
operator norm. In several of the proofs in this 
dissertation it is expeditious to use some specific 
norm, say ||^, for the finite dimensional space involved 
In each such case norm equivalence in finite dimensional 
spaces is envoked to obtain the desired result, i.e., 
there exist numbers , N 2 > 0 such that 

LLl i II i» 2 H r 

For a square matrix A define y [A] by 

r . -1 n . I.I+hAl-1 

y [A] = Urn t- 1 — - , 

h-H-0 

where I is the identity matrix. For a proof that this 
limit always exists and for further discussion of the 
function y see [8], 

For e > 0 and m a positive Integer, let B(e,m) 
represent the open ball of radius e in R m . Let UCR k , 

D 

VCR and suppose that f:UxV-KR . Further suppose that 
all m th order derivatives of f with respect to its 
first k variables exist and are continuous on Uxv and 
that all n th order derivatives of f with respect to its 
last l variables exist and are continuous on Uxv. Then 
f is said to be C m (U) flc n (V) :R^. 
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Let W be any open set in R n and J be an Interval in 
R. Suppose that f, g are C(J) /1 C (W) :B n . Then solutions 
of the initial value problem 

(6) x = f (t ,x) , x(x) = c, 

(7) y = f ( t , y ) + g(t ,y ) , y (x) - (c) , 

are locally unique when (x 5 C) is in J xW. Further if 
x(t) = x(t,x,c) is the solution of (6) then it is well 
known [7] that the matrix $ (t , x, ? ) - x,£) is a 

fundamental solution of the variational equation 

z = f x (t ,x(t ,x, O )z , 

9f . 

where f is the matrix with i.jth element It is 

X a X j 

also known that ■|^-(t , x , ?) = -$(t,x s c)f(t,?), [171- Mow 

oT 

suppose that y (t ) = y (t , x , £ ) is the solution of (7) . 

Then one may write 

|jx(t ,s,y (s) ) = <&(t,s,y(s))[-f(s,y(s))+y(s)] 

= 3>(t ,s,y (s) )g(s,y (s) ) , 


and therefore 
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T 

x(t) = y(t ) + j*(t,s,y(s) )g(s,y(s))ds 
t 

This relation, which is a generalization of the variation 
of constants formula, is due to V. M. Alekseev [1], [2] . 
For further discussion of this relation see [5], [63, [73 
and [193. 



THE UNPERTURBED SYSTEM 


In this section sufficient conditions are given for 
existence of positive and nagative manifolds for the 
unperturbed system. Let £ > 0, m > 0, n > 0 be integers. 
Consider the system 

0 ■ d, 

(8) x = f (t ,x) , 

y = g(t,y), 

where d is a constant ^-vector and f and g satisfy 


(ID 

f is C(R) n C 1 (B(e,m)) :R m for some 

e > 0 ; 

(Hi) 

f (t ,0) = 0, t in R; 


(Iiii) 

there exist a, a > 0 such that 



p[f v (t,x)] < -a|x|“, t in R, x in 

B ( e , m ) ; 

(Hi) 

g is C(R)/1 C 1 (B(ri,n) ) :R n for some 

• f . 

o 

A 

CT 

( II ii ) 

g(t , 0 ) = 0, t in R; and. 


(Iliii) 

there exist B, b > 0 such that 



b | y I 6 < -y[-gy (t ,y ) ] , t in R, y in B (n,n). 


Theorem 1 . Let (Ii) through (Iliii) hold. Then for 
x a fixed number 

S + = { (t ,0,x,y) | 0 arbitrary, | x | <e, y=0, t>r}. 
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is a positive manifold for system (8) and 

S" = { (t ,0 ,x,y ) | 0 arbitrary, x=0, | y | <ri , t£x} , 


is a negative manifold for system (8). Further if 

(0 (t ) ,x (t ) ,y (t ) ) = (0 (t , x , v) ,x(t , x , £ ) , 0 ) is any solution 

of system (8) with (x ,v, C,0) in S + then 


( 9 ) 


x(t,x,c)| < 


sL 


(l+A X, 




for some A > 0 and all t > x; and if ( 0 (t ) ,x (t ) ,y (t ) ) 
= ( 0 (t , x , v ) , 0 ,y (t , x , o ) ) is any solution of (8) with 
(x,v,0,a) in S~ then 


( 10 ) 


y (t ,x,a) | £ 


(l-B|a| B (t-x)) 1/e ’ 


for some B > 0 and all t £ x. 

Proof . It is sufficient to show that solutions of 
x = f(t,x) exist for t _> x and are bounded as in (9), 
and that solutions of y = (g(t,y)) exist ^for t £ x and 
are bounded as in (10). 

From (li) it follows that for each x in R and 
| c | < £ there exists a unique solution x(t) = x(t,x,C) 
with x(x,x,£) = J of X = f(t,x) in a neighborhood V of 
x. So for each t in V and |c| < e one may write 



D R | x(t ) 


lim 

h^+0 


V\ 



= lim 
h++0 



where D R jx(t)| denotes the right derivative of | x | at 
t. Since f (t , 0 ) = 0, t in R, the mean value theorem 
gives 

f (t ,x(t ) ) = f (t ,sx(t ) )ds • x ( t ) . 

D 

Therefore 

Dp | x(t ) 1 £ lim 1 1+h f (t ,sx(t ) )ds I -1) I x(t ) I 
h-*+0 n J x 

| I+hf (t , sx (t ) ) | -1 

<_ lim — ^ — r — — ds ] x (t ) | . 

h->+0 J n 

Let (h, be a sequence of positive numbers with limit 
zero and let 

I I+h . f (t , asx (t ) ) |-1 

J As) = - 1 ' x h • • — , i = 1,2,,.,, s in [0,1]. 

i 

Since 

10,(3)! < ^i+f x (t,sx(t)) -fp|i| 
i i 
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^+f x (t ,sx(t ) )-g— -I 
i 1 


= | f (t ,sx(t ) ) 

X 


it follows that {J i (s)> 1=1 is a uniformly bounded sequence 
sequence of uniformly continuous functions on [0,1] . 
Further lim J, (s) = y [f (t , sx (t ) ) ] . Thus by the 

• 1 X 

X->-oo 

Lebesgue Dominated Convergence Theorem 


r I+hf (t,sx(t)) 1-1 


lim 

n++0 


h 


ds - 


y[f v (t ,sx(t ) ) ]ds , 

X 


and 


D R J x (t)| < 


y[f (t ,sx(t ) )]ds )x(t ) J <_ -a 

X 


sx (t ) | a ds | x ( t ) 


0 


0 


ot+1 


x(t ) 


a+1 


The initial value problem u 
has solution 


au 


a+1 


a+1 


, u ( x ) = | C 


(l+A]c| a (t-T)) 1/a ’ 

where A = Using a standard comparison theorem 

([14], Theorem 4.1, p, 26) one finds that 


x(t ) 


< 


, JSi 

(i+a| c! a (t-x)) 1/a 


t > T . 



Hypothesis (III) implies that for each t in R and 
| a | < n there exists a unique solution y (t ) = y (t , t , a ) 
of y = g(t ,y ) in a neighborhood V of x. So for t In V 
and | ct | < n it follows that 


DjJy (t ) 


_ llrr; I y (t )-hy (t ) 1 — [ y (t ) 
h->+0 h 


> -lim E 
h++0 n 


|(l-h 


g (t , sy (t ) )ds)y (t 


- y ( t ) 


-lim 

h-^+0 


I-hg y (t 


sy(t)) |-1 


ds | y ( t ) 


0 


where D^y(t) denotes the left derivative of |y| at t. 

As above one may apply the Lebesgue Dominated Convergence 
Theorem to show that 


D t 1 y ( t ) 



(t ,sy(t ) ) ]ds | y ( t ) 


sy (t ) 


B ds|y(t)| 


Thus 


D L |y(t)| > ||j|y(t) | B+1 , 


b 6 

and Inequality (10) follows as above with B = 



THE PERTURBED SYSTEM 


Again let SL >_ 0, m > 0 , n > 0 be integers. Consider 
the real system of ordinary differential equations 

0 = d + 0(t , 0 , x,y ) , 

(11) x = f (t ,x) + F(t , 8 ,x,y ) , 

y * g(t ,y) + G(t,8,x,y), 

where d is a constant £- vector and the following 
assumptions hold: 

(Illi) f is C(R)/1 C 1 (B(e 0 ,m) ) :R m for some e Q in (0,1); 
( Illii ) f(t,0) = 0, t in R; 

(Illiii ) there exist a, a > 0 such that 

y[f x (t ,x) ] £ -a|x| a , t in R, x in B( e Q ,m) ; 

(IVi) g is C(R)n (B(n ,n ) ) : R n for some n in (0,1): 

( IVii ) g ( t , 0 ) = 0, g (t,0) = 0, t in R; 

v 

(IViii) there exist $, b > 0 such that 

b|y|® <_ -y[-g (t,y)], t in R, y in B(n,n); 

tj 

(IViv) there exist y, c > 0 such that 

8 2 g. 

max a v ~ a ~ v ~' £ c ly| Y > t in R, y in 

i,j,k dy j dy k 


B(n,n) ; 
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(Vi) 


(Vii) 

(Viii) 


(Viv) 


0, F and G are defined and 

C(R)/1 C 1 (R 2 'xB(e 0 ,m)xB(n,n)):R lc for k = % 9 m 
and n, respectively; 

0, F and G are multiply periodic in 0 with 
period w = (u^, . . . ,w^) , > 0, i = 1,2,..., 

0(t, 6,0,0), F(t, 0,0,0), G(t, 0,0,0) vanish 
for (t , 0 ) in R 1 " 1 "^; and. 


there exist positive numbers p, q , r and L 


such that | 0g | 

• 18,1 

• ! 0 yl 

< L ( 

1 X 

1 + 

|y 1 

) p . 

Fel 

• l p x' 

• l p yl 

< L( 

lx 

1+ 

1 y 1 

I) q , 

l Q 0 1 

• l G x' 

1* ' a yl 

< L( 

lx 

1+ 

|y 1 

D r , 

for 

|x| < 

e o> ly 

1 < h 

• 





£; 


Theorem 2 . Let (Illi) through (Viv) hold with 
y + 1 > a, min{p,q,r } > a. Then for all e sufficiently 
small and r any real number there exists a function v in 
C ( [t , °°) xR^xB ( e ,m) ) : R n , v(t , 0 ,x) multiply periodic in 0 

14-0 

with period w, v(t , 0 , 0 ) = 0 for (t,0) in R , such that 

S + = { (t , 0,x,y ) | 0 arbitrary , | x | <e, y=v(t , 0 ,x) , t_>x} 

£ 

is a positive manifold for system (11). In case i = 0 
min(p ,q,r } is replaced by min{q,r } . 
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This theorem Is the main result of this thesis. The 
proof of Theorem 2 will be given after several results 
are established. First an outline of the proof will be 
presented . 

Let 0 < 6 < min{l,n/e} and for x in R, 0 < e < e Q , 
define 

fl(e) = (v in C([x,«) ) f\ C 1 (R £ xb( e,m) ) :R n | 

v(t ,6 ,x) satisfies (Vli)-(VIiii) } . 

(Vli ) v has multiple period u in 0; 

(Vlii) v(t,0,O) = 0, 9 in R l , t < t; and, 

(Vliii) max sup sup . sup { |v(t,0,x) | , |v fi (t.0,x) | , 

0 in R |x|<e 

t 

I v ( t , 0 , x ) | } < 6 . 

•A 

For v in Q(e) define 0 v (t,0,x) = 0 (t , 0 , x , v (t , 0, x ) ) 
and define F v and G v similarly. Let 
(\p v (t , x , 0 , x ) , £ v (t , x , 0 , x) ) denote solutions of 

0 = d + 0 v (t , 9,x) , 

(12v) . v 

x = f (t,x) + F (t , 0 ,x) , 

with initial condition (0,x) at t = x. When the initial 
conditions are clearly understood the solutions will 
be abbreviated (\|; V (t ) , £ v (t ) ) . It will be shown that an 
operator T, given by 
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0 
» 

$ (t ,s+t ,v(s+t ,^ v (s+t ,t ,0 ,x) ,^ V (s+t ,t ,0 ,x) ) ) 

I 

JO 

xG v (s+t ,^ v (s+t ,t , 0 ,x ) , £ v (s+t ,t , 0 ,x ) )ds , 

3 

where $(t,T,a) = y (t , t , a) and y(t,T,d) is the solution 
of y - g(t,y) with y(T,T,a) = a, is defined on Q(t) into 
fi(e) . 


( 13 ) 

(Tv) ( t,0 } x) = 


Note that 

(14 ) (i|j v (t , x , 0+o),x) , £ v (t ,x,0+w,x) ) 

= (lp V (t ,T,d,x) ,£ V (t,T,6 } x)) + ( 0) , 0 ) , 

implies that (Tv) (t ,0+oj,x) = (Tv) (t ,0 ,x) for v in ft(e) . 

The set n(e) may be thought of as a subset of the 

Banach space of bounded continuous functions, v(t,0,x), 

on V(e) = [ t ,<»)xR^xB(e,m) with multiple period w in 9 and 

norm | |v| | - sup| v(t,0,x)|. It will be shown that for e 
V(e ) 

sufficiently small T is a contraction map on n ( e ) and 

u in 

that there exists a funct ion v e ) such that the operator 
T may be extended to u and Tu = u, The function 
w(t) = u(t , if> v (t ) , £ v (t ) ) can be shown to satisfy the 
differential equation 
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y = g(t ,y ) + G ( t ,\p V (t ) , £ v ( t ) , y ) , 


and it follows that u is the function needed t,o describe 

the positive manifold S + of Theorem 2. 

£ 

The first of the lemmas which will be needed is 
Lemma 3 • Let the hypotheses of Theorem 2 hold. For 
e > 0 sufficiently small, all 0, |x| < e , v in Q(e), 
the solution (i|; v (t ) , £ v (t ) ) of (12v) with initial 
condition (0 ,x) at t * x exists for t _> t and satisfies 


? V (t)| 


< 


(l+A|x| a (t-T)) 1/a ’ 


I^(t)| < M, 

^g(t ) ! , |^(t ) | < M, 


for some positive constants A and M independent of the 
choice of v in fl(e). 

Proof. Let 




1 I 

: £ 0 5 

a 

q-a 

> 

(a+l ) 2 q+1 L 



Then for Jxj < e^ and all 0 solutions (^ v (t ) , ? v (t ) ) of 
(12v) exist and satisfy | £ v (t ) | < £j in a neighborhood U 
of t. So for t in U it follows that 
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lim 

h~>+0 


lim 

h-H-0 


r V iV c v 
h 

U v + h(f(t,g v )4P v (t^ v ^ v ))|-lg V 

h 


< lim lkl£l + ]F v (t^ Y ,5 v )|. 

h-H-0 


As in the proof of Theorem 1 one finds that 


lim 

h^+Q 


|g v +hf(t,g v )l-l g v | 
h 




a+1 


From (Vlil) and .(Villi) it follows that 


|v(t,e,3C) ] < 5 | X 1 £ | X I , 

and from (Viii) and (Viv) 

lp v (t,* v J e v >] < 2 l?+1 L|5 T l q+1 . 


Therefore 

d r I I sif-i s v 1 01+1 + 2 q+1 L|sY q+1 

<_^( a - 2 1«(a + X)L £ q - a )|5 V | a+1 

i=||E V |“ +1 


5 



where A 


2 

= ^(a-(a+l)2 q+1 Le q_a ) > 0. Thus, as in 
Theorem 1, 


£ v (t ,x,9,x) | 


( 1+.A |x| a (t-x)) 1/a ’ 


T . 


The existence and continuity of £ 0 , and 

in V(e, ) are well known, see ([17] a p. 25). Let 

A X 

( X ( t ) , X ( t ) ) = Og-(t),cJ(t)) or C^ 1 x ( t ) ,?x(t ) ) , then 
from (12) it can be seen that 


X = O 0 +e y v e )x + (Q x +0 y v x )x, 
x - (f x + F x +Vx )x + (P 0 +? y v 9 )Xj 


where 0 Q - 0 Q (t , t|/ v (t , x , 9 ,x ) , £ v (t , x , 0 ,x ) , v (t , (t , t , 6 , x ) , 
£ V (t,x,9,x) ) ) , 
etc. Therefore 

| X I <_ L[ ( | C V | + |v (t , 'p v ,C V ) | ) p + ( | ? V | + j v (t , if) v , £ v ) I ) P 

iv 6 mxi 

+ L[ ( ) £ V j + J v (t , \p V , £ V ) | ) P + ( j C V I + I v (t , \p V , £ V ) | ) p 
I v x I ] I X I 

< 2 P+1 L | ,§ V j P | X ( + 2 P+1 LU V |P]x |, 


and D R jxj = lim 1 X+hx I - i X 

h-H-0 h 5 


\X* 



24 

| I+hf |-1 

< lim 7 — Ixl + Ip +F v j |xl + |f q +f v a I I X | 

- h _^ +0 h 1 A 1 1 X y x 1 1 A 1 1 6 y 0 M 1 

i yCf x (t,g v )]|xl + 2* +x L\e\*\x\ + 2^ 1 LU v | q Ul 

< -A 0 U V | a |x| +2 q+1 L|C V | q |x|, 

where A Q = a - 2 q+1 Le q " a > 0. Since D R | X | < | A | and 

00 

since p > a implies |£ v (s)| p ds < 00 it follows that there 

* 

T 

exist constants >0 such that 

t 

(15) | X ( t ) | < K 1 + K 2 jlx(s) | |C V (s) | p ds, 

T 

so that 

D r | X ( t ) | < -A 0 |C V (t)| a | X (t)| + 2 q+1 L|5 V | q 

t 

x(K 1 +K 2 I X ( s ) I U V (s) | P ds ) . 

t 

Now consider the scalar equation 

v = -A 0 |C V | a v + K 3 U V | q , 

where 

op 

0 < K 3 = 2 q+1 L(K 1 +K 2 2|? v (s)| P ds) 

T 


< *>. 
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It is clear that for t _> x 

t 

-v i a. 


v(t) = v(x)expC-A 0 j | ^ |°as] + | K | q exp[-A Q 

t T T 

I 5 V (s) | q ds. 


U V | a du]ds, 


< V ( T ) + K- 


Now for 


e < 


A(q-a) 


q-a 


and ( x | < e. 


K 


| £ V (s ,x, 9,x) | q ds £ K, 


* 3 ° 


1 x I q ds 


K 3 a|x| q - a 


(1+A |x | a (s-T) ) q//a A(q-a) 

l l 

Then if 0 < v(x) £ 1 it follows that 0 £ v(x) £ 2 for 
all t £ t. Suppose v(x) is chosen so that v(x) = j X ( ^ ) i 
then |x(t)| £ v(t) for all t _> x. For suppose this is 
not the case. Then there exist t £ x such that 
| X (t ) | > v(t ) . Let 


t, = sup {t | | x(t ) | — V ( t )£0 , X£t £t } . 

-1* A 

t 


then j x Ct ) | = v (t x ) and 
D R lx(t l )| £ -A 0 | 5 v (t 1 )|%(t 1 ) 

t 

+ 2 q+1 L|C V (t 1 ) | q (K 1 +K 2 .[ U V (s) | P v(s)ds) 

4 

-A 0 |C V (t 1 ) | a v(t 1 ) + K 3 |C V (t 1 )| q = D R v(t 1 ) 
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Thus 

D R (| X (t 1 )|-v(t 1 )) < 0, 

and this contradicts the definition of t^. Therefore 

from (15 ) it follows that 

t 

| X (t ) | < K + K 2 2 |5 V (s) | P ds 

T 

00 

£ K 1 + 2K 2 | u V (s) | P ds = M < «° . 

T 

Lemma 4 . Let ( I Vi ) through (IViv) hold and let 
y(t,x ,a) be the solution of y = g(t,y) with y (x, x,a) = a. 
Denote -^-yCt^a) by $(t ,x,a) . Then for |a| < n and 
t £ x 

| $ (t , x, a ) | < 1, 

and 

|^-$(t ,t,o) | £ T | t-x | | a | T > 0, j =1,2,... . 

3 

Proof . For solutions z(t) of the linear system 

z = C ( t)z , 

it is known that ([18], Theorem 3, p. 58) for t £ x 

t 

| z (t ) | £ | z ( x ) | exp j— y[-C(s)]ds. 

x 

Since $(t,x,c) is a fundamental matrix for 
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z = g y (t,y(t,x,a))z, 
it follows that 

t 

j $(t , x,a) | < exp -y[-g (s,y(s,x,a))]ds 

x y 

T 

<_ exp -b | y ( s , t , a ) j ^ds <_ 1 . 

« 

t 

Since g is C (R) /I C 2 (B(n,n) ) :R n the partial 
derivatives of $(t ,t,cf) with respect to each component 

c. of cr exist and satisfy 

J 

♦ 

Z = g y (t ,y(t ,x,ff))Z + A. (t ) , 

where 

A^(t) = 1YL (t )$ (t ,x,a) , 

and Mj(t) is the n x n matrix with k, £ th element given 
by 

n 3 2 g, 3y, 

Ji8775y^ !t (t )T • 3 > ’ T ’ a } * 

See ( [14 ] , Theorem 3 • 1 » P* 95). Thus , by the usual 
procedure for solving a nonhomogene ous system of linear 
differential equations ([14], pp. 48—49) it can be seen 
that 

t 

|4~ (t,x,a) = $ ( t , s , o ) A . ( s ) d s . 

j J J 


T 



Therefore for t < r 
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3 a 


j 


-4> (t , tj a) 


$(t ,s ,a) 


|A j( s) 


ds < 


A . (s ) ds 
3 


Making use of norm equivalence in R n there exists N > 0 
such that 


| A . ( s ) | <_ N max 
9 i,k, £ 


a .2 

9 Si 


9y i 9y £ 
£ Nc | y (s jTjff) | Y . 


(s,y(s,x,a) ) 


f(s,T,a) 


Therefore by inequality (10 ) 


3 


< Nc 


y (s , x , a ) | Y ds £Nc | a | Y 


( 1-B a y (s-t)) y 


7 $ £ 


< Nc I a l Y (T-t ) . 


Choose r * Nc and the result follows . 

Theorem 5 • Let the hypotheses of Theorem 2 hold. Then 
for e > 0 sufficiently small the operator T given by 
equation (13) is defined on 0(e) and maps Q(e) into 
itself. 

Proof s Let > 0 be small enough so that Lemmas 3 and 
A hold. Then for (t ,0 ,x) in V.Ceg) 



oo 


29 


Tv(t ,6 ,x) | < 


L( I c v (t+s ,t ,0 ,x) 


0 


+ | v (t+s ,^ v (t+s ,t ,0 ,x) , £ v (t+s , t ,0 ,.x) ) | ) r ’ rx ds 


,r+l 


< 2 r+1 L 


£ v (t+s , t , 0 ,x ) J r+1 ds 


0 

OO 


< 2 r+1 L 


X 


r+1 


(1+A x a s )(r+ 1 ) / a 


ds 


. 2 rtl toU [ rtl ' 11 
A (r+l-a) 


< OO 


Let 

j. i 

r+l-a I 

r 

then for e < T is defined on fl(e), Tv (t , 0,0) = 0 
and | Tv (t , 0 ,x ) | < 6. Also since 


= min < e 2 , 


A6 (r+l-a) 


2 r+1 La 


x 


i r+1 


(1+A x 


a g ^ (r+1) /a 


( 1+As ) 


(r+1 )/a 


and 


00 

V" ■ , v ) ds < 00 It follows that Tv converges 

J(l+As) (r+1}/a 


uniformly for (t , 0,x) In V(e) which Implies that Tv 

is continuous on V( e ) . It has been noted previously 

that Tv has period w in 0 for v in ft(e). It remains to 
3 3 

be shown that ^tt(Tv) — -(Tv) exist and are continuous 

d b d X 

and bounded by 6 for (t,0,x) in V( e) . 
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In the following the summation convention will be 
used. Let A represent either the 0 or the x vector. 
Then the i,jth element of the matrix - 5 — ( $G ) is given 

o X 

by 


- — ($., gJ) 

3 X, ik k' 


3$ 


ik 


3o 


l 

+ $ 


3 V £ 3 $. 
38 


y v 

m + 


G 


ik 

3^ 


sy 


m 9X j 

3* 

m 9X 

9 J 

3G k 3^ 

4 . 

9G k 

35 

"a 3 V 


3x * 

3y 

' 9V £ 

4* 

8 h 

h; 

30 m 3Aj 


3x 

m 

3 A 


where 


$G V 


$(t ,t+S , v(t+s , ^ v (t+s ,t , 0 ,x) , £ V (t+S ,t , 0 ,x) ) ) 
xG v (t+s ,i(j V (t+S ,t , 0 ,x) , £ V (t+s,t , 0 ,x) ) . 


From Lemmas 3 and 4, (Viii), (Viv) and norm equivalence 
for R n there exists N* > 0 such that 


-|yUG V D < NN*r(s) Jv| Y (2M)L( |g v |+jv| ) r+1 

+ N*[L( | ^ | + I v I ) r M+L ( | 5 V I + j v | ) r 3M ] 

s|x|Y +r+1 

£ Mo (l+A|x| a s) ( Y +r+1)/ a 

+ M - 1-1 - • ■ , 

0 (1+A | x | a s ) r/a 


where Mq = 2 r N*LM max{4Nr,3). Since r > max{2a-(y+l) ,a} 


if follows that 
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s x 


y+r+1 


rdLs = 


. . 2.1 . |.y+r.+l-2.a 
a x ' 


0 


(1+A x s) 


a. N ( y+r+r J7a «2 


A (a- (y+r+1) ) (2 a- (y+r+1) ) 


x 


x 


ir-a 


(1+A x s ) 


a sr/a ~ A (r+l-a) 3 


and 


oo 


(l + A]x|“s) (Y+r+1)/a ^ S " |(l+As) ^ +r+1)/a_ ^ S < 


t y+r+X 


x 


0 


7 — ds < 

(1+A lx “s) r/a - 


0 


( 1 + A s ) 


r/a 


-ds < 00 . 


Q 


Therefore ■| T (4>G v )ds converges uniformly for (t , 9 , x ) in 
j d A 

OO 

V ( Eg ) and is dominated by a continuous real valued 
function D(e) defined on [0, £p ) with lim D(e) = 0. If 


e-»-0 


> 0 is such that D(e^) < 5 then for (t,8,x) in V( ) 


1T ( Tv > - ll (TCV)ds > 


and 


fx (T '' ) 


< 6. 


The following lemmas will be used to show that T 
is a contraction map on Q(e), 
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Lemma 6 . Let the hypotheses of Theorem 2 hold and let 
> 0 be such that Theorem 5 holds. Further let 
v, w be in n(e^) and suppose that 

(^ V (t , T , 0 s x) ,£ V (t ,T, 0 , x ) ) , (> V (t , T , 0 , X ) , (t , T , 0 , X ) ) 

are the corresponding solutions of (12v), ( 12w ) . Then 
there exist K > 0, independent of v, w in such 

that 


U V -^ W | + |C V -? W | < K | |v-w| I, 


where j j v-w | | = sup | v ( t , t , x ) -w ( t , 0 , x ) | . 

V(e 4 ) 

AAA 

Proof . For (t,0,x,y) , (t,6,x,y) in V(e^ ) xB (n ,n ) and 


x | = max { | x | , | x | } , |y| = max { | y | 3 | y | } , 


it is clear that 


I 0 ( t j 0 , x , y ) - 0 ( t , 0 , x , y ) | £ L( |x | + |y | ) p ( | 6-0 | + |x-x | + |y-y | ) , 
Let h(t) = max{ | £ v (t ) | , | £ w (t ) | }, then 

max{ | v(t , ij; v (t ) , ? v (t ) ) | , |w(t ,^ w (t ) , C w (t ) ) | } <_ h(t ) . But 
| v(t ,4) v ^ v )-w(t | 

< | V ( t j i/j v , € v ) -v ( t , ) | + I v ( t , ) ~w ( t , ) | 

< ( |^ V -*|» W |+U V -£ W |+| | v-w | | ), 
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implying that 

1 0 V (t ,if> v , ? v )-0 w (t , £ w ) j 

< 2 p Lh p (t)(2( |^ V -^ W | + U V -C W | )+| | v— w | |) 

< 2 p+1 Lh p ( | ip v -ip W | + | £ V -£ W | + | | v-w | | ) . 

In the same way, 

|P v (t,^ v ,5 v )-P w (t,^,5 w ) | 

< 2 q+1 Lh q ( t ) ( 1 * V -* W I + U V -S W 1 + I i v-w j I ) . 

Let H(t ) = Lh(t) max{2 p+1 ,2 q+1 >. Note that Lemma 3 

OO CO 

implies H p (s)ds < 00 and H q (s)ds < 

J J 

T x 

Now from (12) 

(16) D r |/-^ w | < | F (^ v -i); w ) 

< |0 v (t,i|; v ,C V )-0 W (t^ w ,? w )| 

< H p ( t) ( U V -^ W | + U V -C W | + ! | V-w I i). 


and 
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-V ,W| 


D |e v -£ W ! = lim J .. ^ -. g .. ) +h ^ -1 U-ll -k 
h-*+0 h 

< lim l(g V -g W )+h(f (t;e V )-f (t ,g W )) klg V - 


h-^+O h 

, I mV/, , V r V N „W / , ,W ,_W A I 

+ I F (t,| ,?■ )-P (t ,1{J ,£ ) I . 


The mean value theorem gives, for some v in [ 0 , 


f(t,? v ) - f(t,5 w ) = f x (t , v£ v + ( 1-v ) 5 W ) ( £ v -£ w ) , 


and therefore 


D R |C V -5 W | < y[f x (t,v5 v +(l-v)C w )]|5 v -C , 'l 

+ H q (t)(|i|) V -* w | + |e T -c w | + | |v-w| |). 


But 


u[f x (t ,v? v +(l-v)? w )] < -a|v? v +(l-v)5 w | a < 0, 
so that 


d rU V "^ W | £ H q (t)( 1/-^ W | + U V -5 W | + 1 )v-w| |), t 
The above inequality, together with (16), gives 
d r (|^ w | + |s v -s w |) 1 K(t)(|^ w | + U v -C w l + I |v 



> T , 


-w | | ) , 
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where K(t) = (t ) + H^(t). Since j ^ ( t)-^ W (t) | = 0 
= | £ v ( t)-£ w ( t) I it follows that 


,v ,w . r v r w 
ip -ip + 5 -C £ v-w 


K(s) exp { 


K(u)du }ds 


K | | v-w ] | , 


where K 


exp 


K(s )ds 


< 00 . 


X 


Lemma 7 • Let the hypotheses of Theorem 2 hold and let 
$(t , x,o) be as in Lemma 4. Then for all (t , 0 , x ,y ) , 

AAA n 

(t,0,x,y) in (- 00 , x ] xR xB ( ,m) xB( n ,m) there exists P > 0 

such that 


AAA 


l$(t 3 x J y)G(t,e j x 3 y)-$(t , x,y )G (t 3 0 3 x 3 y) | 

< P{( |x |+ |y | ) r +(x-t ) ( |x |+ |y | ) r+1 ! y | Y } 

A A 'N 

x{|0-0|+|x-x|+|y-y | }, 

/N A 

where | x | = max { | x 1 3 j x | } , | y | = max { | y | , ] y | } • 

Proof . The mean value theorem gives 


AAA 


$(t 3 x 3 y )G(t 3 0 3 x 3 y )-$(t 3 x 3 y)(Ut 3 0 3 x 3 y ) 


sup 




9x 


( $G) 


3y 


($G) 


{ | 0-0 | + | x-x | + | y-y | } 3 


where the supremum is taken over R xb ( e^ 3 m)xB(n,m) 
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For X = 9 or x it follows from Lemma 4 and (Viv) that 

|y($ (t 3 t , y )G(t ,Q ,x,y ) ) < L( | x | + |y ] ) r . 

n 

Again norm equivalence for R implies there exists > 0 
such that 

By applying Lemma 4 and (Viv) it can be seen that 

|y($(t 3 T,y)G(t,0,x,y) < N^L( |x|+|y | ) r+1 r | y | Y ( t - t ) 

+ L( | xj+jy | ) r . 

Put P = max{N^Lr,3L} and the result follows. 

Theorem 8 . Let the hypotheses of Theorem 2 hold. Then 

for all e. > 0 sufficiently small the operator T, defined 

by equation (13) > is a contraction map on fl ( e ) . 

Proof. Let e c > 0 be such that Theorem 5 holds and let 
. b 

v, w be in Q(e^). Then for h(t ) = max{ | £ v (t ) | , | £ w (t ) | } 
it follows from Lemmas 3 , 6 and 7 and equation (13) 


j t , y )G (t , 0 jX,y ) 

< N x | G | max |—(t ,x,y) + | $ (t , x ,y ) | ~G (t , 0 , x ,y ) 

J 2 


that 
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Tv-Tw < 


$ (t ,t+s , v(t+s ,^ v (t+s ,t , 0 jX) , £ v (t+s a t , 0 ,x ) ) ) 


*G V (t+s ,^ V (t+s ,t , 9 ,x) , £ v (t+s ,t , 9 ,x) ) 

- $ (t ,t+s ,w(t+s ,^ W (t+s ,t , 6 ,x) , £ w (t+s ,t , 9 ,x) ) ) 
*G w (t +s ,\}> W (t +s ,t , 6 , x) , £ w (t+s ,t , 0 jX) ) ) | as 


< P j C 2 r h r (s )+2 r+1 sh r+ ^ +1 (s ) ) 


x ( j ) + j £ v -£ w | + | j v-w | | )ds 


< 2 r+1 P(K+l) 


(h r (s )+sh r+Y+1 (s ) )d 


s v-w 


< 2 r+1 P(K+l) 


00 

r 


S X 


(l+A!x| a s) r/a 
l r 

ds 


ds 


0 

,r+l 


(l+A|x| a s) Tr+ ‘ Y+l)/ 


a 


v-w 


< 2 P ( K+l )D ( I x ) j v-w , 


where D ( e ) is described in Theorem 5- Let Eg > 0 be 
such that 2 r+1 P (K+l )D ( Eg ) < 1. Then for e = mint£g,£g} 
the result follows. 


Proof of Theorem 2 . Since T maps 9,(e) into itself and 
is a contraction map it follows that the sequence 

{v n } n=l’ V 1 in v 2 = Tv l j *“’ v n+1 = Tv n 5 * * * ls 

Cauchy and therefore converges to a continuous function 

u in fi( e ) . Since {v n } n °*^ uniformly Lipschitz with 

Lipschitz constant 6 it follows that the limit function 



36 


u is likewise Lipsehitz. By extending: the definitions 
of 0 V , G v and F v to allow v to be continuous it can be 
seen that the initial value problem (12u) has the 
unique solution <> u (t , t , 0 ,x ) , £ u (.t , t , 0 ,x ) ) . Further, 

V . v 

since 0 n (t ,9 ,x) 0 u (t,j,x) and F n (t,6,x) F u U,8,x; 

uniformly on [T,~)xR K 'xB(e } m)xB(n,n) it follows 

v v 

([14], Theorem 2.4, p. 4) that (i p n ,% n ) •> (^ u ,^ u ) 

a 

uniformly on IxJxR xB(e,m) where I is any compact 

interval in Observe that the bounds computed 

for £ v in Lemma 3 hold for £ u . 

Consider the sequence of functions (J } °°, , where 

n n=l’ 

v v 

J n (s) = f(t,t+s,v (t+s,ij> (t+s,t,e,x) ,c n (t+s,t,9,x) )) 

V V V 

xg n (t+s,^ n (t+s,t,e,x),£ n (t+s,t ,e,x) ) , 

n = 1,2,.. . 

From the above remarks and from the continuity of 
and G it follows that 

lim J ( s ) = $(t,t+s,u(t+s,i U (t+s,t ,9,x),C U (t+s,t,6,x) )] 

n+oo 

xG U (t+s ,^ u (t+s ,t ,9 ,x) , £ u (t+s,t ,9 ,x) ) , 
f o r x < t < s < 00 . 
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It follows from Lemmas 3, 6 and (Vlv) that 


v 


v 


v 


J n (s) I < j G n ( t +s , ^ n ( t +s , t , 0 , X ) , 5 n (t+s,t ,0,x) ) I , 


v 0 r+l T i ir+1 

£ 2 r+x L|C n (t+s,t,e,x) | < 2 L ' x ' 


(1+Al.x a s) r+1/a5 


n = 1 , 2 ,...,, 


therefore the Lebesgue Dominated Convergence Theorem 
implies 


(18) u{t , 0 ,x ) 


lim v ( t , 6 , x ) 

n^-oo 


0 


= lim 
n-*- 00 


J n _i Cs )ds 


0 

lim J n _^(s)ds 

i n + ~ n 


0 

$ ( t , t +s s u ( t +s j 4> u ( t +s , t , 0 , x ) , 
£ u (t+s ,t ,0,x) ) ) 


xG u (t+s 3 \p U (t+s ,t , 0 , x) , £ u (t+s , t , 0 ,.x) )ds 


By using this representation of u it will now be shown 
that the function u(t) = u (t ,ij> u (t ) , £ u (t ) ) satisfies the 
differential equation 


y = g(t ,y ) + G(t,^ u (t),C u (t),y). 


From the definition of w(t) it follows that 



0 


(o(t) = 


$(t ,t+S ,u(t+S , l|) U (t+S ,t , ^ U (t ) , C U (t ) ) 3 

£ U (t +s ,t , if> U (t ) , £ U (t ) ) ) 
xG u (t+s ,^ U (t+s ,t ,^ u (t ) , £ U (t ) ) , 
C u (t+s,t ,^ U (t ) ,£ u (t ) )ds , 
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but from the uniqueness of solutions of (12u), 

(5 u (t+s,t ,\|> u Ct > x J 0 3 x),5 u (t,T,e j x)) 3 
C u (t+s,t ,^ u (t ,T,e,x) s £ u (t 5 T,e,x) ) 

= (ljj U (t+S jT^jX) , 5 U (t+S,T,e,x) ) • 


Thus 


0 

U)(t ) = j$ (t ,t+S ,u (t+s , !^ U (t+S , T) 0 ,x) , £ U (t+S , T, 0 ,x) ) ) 

°°xG U (t +S ,^ U (t+S,T,0jX),^ U (t+S,Tj0,x))dS 
t 

= I'l>(t 3 s j u(s,^ u (s),? u (s))) 

°°XG ( s , \jj u ( s ) , £ u ( S ) ,u ( s , T^ u ( s ) £ u ( S ) ) )ds 


$(t,s,w(s))G(s,^ u (s),? u (s),co(s))ds. 


Let t be fixed such that t < t < °° and let 


Q fc (t) 


t 

$(t ,s,w(s) )G(s,^ u (s) ,5 U (s) ,w(s) )ds, 
t 


and 



Q~ (t ) 


f 

$ (t 3 s v o> ( s ) ) G ( s 5 4* u ( s ) , £ u ( s ) , a) (s ) )ds . 
t 
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Since ^r(t 3 s,a) = g y (t >y (t 5 s , a ) ) $ (t , s , a ) is continuous 
for t £ s £ t and ]a| < n and since $ (t ,cr) = I 




-G (t , i^ U (t ) , C U (t ) j w(t ) ) 


g y (t ,y (t ,s,w(s) ) )$(t ,s ,w(s) ) 


*G ( s , tp u ( s ) , £ U ( s ) , co ( s ) )ds 


Also since 


g y (t ,y (t ,s ,w(s ) ) )$(t ,s , to ( s ) ) 


*G ( s , ( s ) , £ u ( s ) , u>( s ) )ds 


|y(t jSjoKs) ) I Y+1 ( |? u (s) | + |w(s) I ) r+1 ds 


2 r+1 Lc 

Y+l 


|? U (s)| r+Y+2 d3 


. "I 00 

2 r+1 Lc ! 


< * x 

- Y+1 ia+A( S -T)) (r,+Y+2:i/a 


7— ds 


implies that 





00 


g y (t ,y(t ,s ,w(s) ) )$(t ,s ,w(.p) ) 
xG (s , (s ) , £ u (s ) , go(s ) )ds , 


It follows that 


(20) w(t) = 


at Qt(t) + lt Q t (t) 


= -G(t J ^ U (t),C U (t),o)(t)) 


g y (t » y (t , s , w( s ) ) ) $ (t , s , w(s ) 


xG (s , ^ u (s ) , £ u (s ) , w(s ) )ds 


u , 


Now observe that 
t 

* J 

jgg(t ,y (t ,s, w(s ) ) )ds = g(t , <o(t ) ) 

00 

- 11m g (t ,y (t ,s,w(s)) . 

S-VOO 

S >t 

But 


11m |y(t,s,w(s)) | < 11m — = o, 

3+°° s-*°° (l-B(t-s) | «(s) | p ) 1/p 

s >t s >t 


and since g(t,0) = 0, t In R, and g is continuous. It 
follows that 11m g(t ,y (t ,s,io(s) ) = 0. Also 

S-^oo 

s>t 
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(t ,y (t jS 5 (a)(s ) ) ) — Sy ( t >y(t ,y (t ,s ,co(s ) ) ) 

■|~(t 9 s y o)(s ) ) >s } w(s) ) w(s ) 

- g (t ,y ( t , s , w( s ) ) « (t 5 s , u)( s ) ) 
j 

x[-g(s,w(s) )+<I>(s)]. 

Therefore 

t 

g(t j w(t ) ) = g (t ,y(t ,s,w(s)))$(t ,s,w(s) 

J J 
oo 

x [ cL ( s )-g(s , w(s ) ) Ids . 

This equation, together with (20), gives 
a) ( t ) - g(t ,co(t ) ) - G (t , \p u (t ) , £ U (t ) , u)(t ) ) 

OO 

= - g (t,y(t,s,w(s)))f(t,s,u(s)) 

j t y 

t 

x[ oo( s )-g(s , w(s ) )-G (s , (s ) , £ U (s ) , to(s ) ) ]as 

Let 4>(t ) = d)(t) - g(t , to ( t ) ) - G(t ,i|> u (t ) ,£ U (t ) , a(t ) ) , 
then <j) satisfies 

OO 

(21) 4> (t ) = g (t,y(t,s,w(s))$(t ,s,w(s))<j>(s)ds. 

v 

t 
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It- will be shown that <J) (t ) = 0, t in [x,°°), is the onlj, 
solution of (21) and the proof of Theorem 2 will be 
complete. First observe that <j) (t ) is uniformly bounded 
on [x,°°). For it has already been shown that 
u)(t) + G(t ,^ U (t ) ,£ u (t ) ,u)(t ) ) Is uniformly bounded, see 
(19) and (20), and (IVii), (IViv) and the definition 
of w imply that g(t,w(t)) Is uniformly bounded. 

From Lemmas 3, (IVii), (IViv) and the definition 

of w 


(j) ( t ) | < 


Y+l 


| w(s ) | Y+l | cf>(s ) |ds £ ^-f |C U (s) | T+1 1 <() ( s ) |ds 


oo 

St{ ( v+ ' l V7ct 1 ^ ( 5 } I d s 

Y+1 i (l+A(s-r)) lY+1)/a 


( Y+l ) ( Y+l-a) (1+A(t _ T) )(Y+l-a)/a sup|4,(sj|. 


Choose t £ x so large that 

(Y+l)(Y+l-cO 7 l+A(t^T)) (Y+1 " a)/a < X ' ^ h6n f ° r a11 

t > t 


| <f> ( t ) | < sup 1 4 > ( s ) | , 

S£t 

but this Implies that <J>(t) = 0, t > t. Thus 



/V 
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|4>(t)| 



| <p (s ) |ds 


5 


and therefore 


d 

dt 



L'| 


(pis) |ds 


< 0 . 


But this gives 


t 

’ /N 

| <(> ( s > | ds >_ 0 , t £ t , 
t 


/s 

and therefore <J>(s) = 0, t < s < t. 


Corollary 9 « Solutions of 

y = g(t ,y ) + G(t ,ijj u (t ) ,£ u (t ) ,y) on S* are 0(t -1/a ) as 
t °°. 

Corollary 10 . The positive manifold S* is unique. 

Proof . Equation ( 18 ) shows that T may be extended to 
9, ( e ) and since Lemma 6 holds in sTCiT the result follows 
from inequality (17). 

The corresponding theorem for negative manifolds 
will be stated here. Consider the following hypotheses. 

(VIII) f is C(R)OC 2 (B(e,m)):R m for some e in 


( 0 , 1 ); 
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(Villi) g ( t , 0 ) = 0, f x (t ,0) = 0, t in R; 

(V'lliii ) there exist a, a > 0 such that 

y[f (t ,s) ] < -a | x | a , t in R, x in B(e,m); 
(VIliv ) there exist tc, K > 0 such that 

3f, 

max (t ,x) <_ K I x | K a t in R, 

i j k 

x in B(e,m) ; 

(Villi ) g is C (R ) A (B(ri,n) ) :R n for some r\ in 

( 0 , 1 ); 

(Vlllii ) g(t ,0) = 0, t in R; and, 

(Vllliii ) there exist 6, b > 0 such that 

b|y| e < -y[-g y (t ,y ) ] , t in R, y in B(n,n). 


Theorem 11 . Let (Vlli) through (Vllliii) and (Vi) 
through (Viv ) hold with K + 1 > 3 and min{p,q,r} > 3. 
Then for n > 0 sufficiently small and t any real 
number there exists a unique function w in 
C( (-«> s T]xR^xB(n,m) ) :R m ,w(t ,0 ,y ) multiply periodic in 
6 with period w, w(t , 0 ,0 ) = 0 for (t ,0) in (- 00 , t]xr^ s 
such that 


S” = { (t ,0,x,y) | 0 arbitrary, x=w(t ,0,y), |y | <n, t<t}, 

is a negative manifold for system (11). If the 
0-equation is absent then min{p ,q,r } = min{q,r } . 



The proof of Theorem 11 follows closely that of 
Theorem 2 with the exception that one considers t < x 
Instead of t > x, For w in the proper class of 
functions one proceeds to prove Theorem 11 by finding 
a solution for the system 


( 2 Iw ) 


0 = d + 0(t,0,w(t,6,y)j y ) , 
y - g(t ,y) + G (t , 6 ,w (t , 0 ,y ) ,y ) , 


w(t ,0,x) 


¥ ( t , t +s , w ( t +s , cj) W ( t +s , t , 0 , y ) v w (t +s , t , 0 , y ) )) 

0 

xF(t+s ,tj) W (t+s ,t ,0,y),v w (t+s,t j 0 ,y ) , 
w(t+s ,cj) W (t+s ,t ,0 ,y ) ,v w (t+s,t , 0 ,y ) ) )ds , 


where ( <j> w (t , x , 0 ,y ) , v w (t , x , 0 ,y ) ) is the unique solut ion 

of (21w) with initial condition (0,y) at t = x and 

¥(t,x,c) = -5— -(t,t,?) } x(tjXjC) is the unique solution 
0 Q 

of x = f(t,x) satisfying the initial condition z, at 
t = x . 


Corollary 12 . Solutions of 

x = f(t,x) + F (t , q w (t ) s x , v w (t ) ) are 0(Jtj“^ / ^) as t 

Remark . The hypothesis of Theorem 2 which restricts 
y + 1 > a can be relaxed in some cases to include 
Y + 1 = a. 
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Suppose that f(t,x) satisfies (Illi) through 

(llliii) with a = y + 1. Let be a sequence of 

positive numbers with < a and lim a, = a. Suppose 

v k-* 00 

there exist s a sequence {f (t ,x) of 

C (R) 0 C 1 (B( e »m) ) :R m functions with f^(t,0) = 0, 

ot 

t in R, and y[f (t ,x) ] < —a | x | k for each 

k = 1 , 2 ,*.. . Finally suppose that lim f^ k; (t,x) 

k-*°° 

= f ( t , x ) uniformly on IxB(e,m) , I any compact subset 
of R. 

By applying Theorem 2 it can be seen that the 
system 

e = d + 0(t ,0,x,y) , 

(22k) x = f (k) (t,x) + F ( t , 0 , x , y ) , 

y = g(t ,y) + G(t,0,x,y), 

has a positive manifold 

S £ | k = { (t ,0 ,x,y ) | 0 arbitrary, |x | <e, y=u k (t , 0 ,x) ,t> r } , 
provided min{p,q,r} > a and e sufficiently small. 
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Let (^(t ) , ? k (t ) ) represent the unique solution 
of 


(23k) 



e = d + 0(t,0,x,u k (t ,0,x)), 
fC^CtjX) + P(t,0,x,u k (t,0,x)), 


with initial value (0 ,x) at t » t. Then 

w k (t ) = u k (t ,^ k (t ) j? k (t ) ) is the unique solution of 

y " g(t,y) + G(t,^ k (t),? k (t ),y). 


with initial value u k (x,0,x) at t = t. 
Since the functions 


© k (t) * d + 0(t ,j^ k (t ) .,S k (t ) ,.a» k (t ) ) , 

P k (t) = f (k) (t,C k (t)) + P(t ,^ k (t),C k (t) J w k (t)), 
G k (t) = g(t,w k (t)) + G(t , 4y(t ) ,C k (t),co k (t)), 

are uniformly bounded on [x ,°°) it follows that the 
sequence { ( ^ ^ 5- k » w k ) } k -i i ;S uniformly Lipschitz and 
therefore equicontinuous . Now for each k * 1,2, , 
( 4> k ( x ) , S k ( X ) , w k ( x ) ) = (0,x,u k (x,0 ,x) ) . But 
{u k ( x, 0 , x ) ) k ” 1 is a bounded sequence and therefore 
has a limit point , say u Q . Let ( u k (^) (t,0,x) be 

a convergent subsequence of {u k (x,0,x) with 

lim U k(1) (x, 9 ,x) = u Q . Since { (\» w k ) > k * x is 



5o 


equieontinuous and since solutions of (22k) are unique 
if follows that { 1^* 5j_> ) ^i=i» where i = k(i;, 

converges uniformly on [t,°°) to a function (^q,£ 0 ,cu~) . 
Therefore, from the continuity of © k , F k and G k it 
follows that { (0 k ,F k ,G k ) > + (0 o ,F o ,G o ) uniformly 
on [t,°°). Thus (i|jQ,5Q,ajQ) is the unique solution of 
(22o) with initial value (0,x,Uq) at t = x. 

Observe that the asymptotic behavior of the 
solutions of the x and y equations in the limit is the 
same as in the case a < y + 1. For it is clear that 

U k (t,t,e,x) | < 

N k Ct)| < 6U k (t) j 

for some A k > 0, k = 1,2,... . Since A k A, a k a 
as k + « it follows that 

U o (t,I,0,x)| < 

| w Q (t ) | <_ 6 | C 0 ( t ) | . 


(1+A x| a (t-T)) 1/a 
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in case system (22k) is autonomous it is clear 
that (u^.(t ,0 ,x) } - (u k (0 ,x) } converges uniformly on 

o 

R xB(e ,m) . Thus for autonomous systems one may obtain 
a representation of the form Uq = TUq for the integral 
manifold even in the case y + 1 = a. 

In view of the above discussion it can be seen 
that Theorem 11 can be extended to include the case 
0 = k + 1. 

Now it can be seen that the model problem dis- 
cussed in the introduction has the predicted solution. 
In that case a = 0 = 2 and y = k - 1* So if X and Y 

q 

are o( ( |x| + |y | )~ > ) one may apply an extended version 

3 _1 

of Theorem 2 by letting f^ k ^(t,x) = -x k to obtain 
a positive manifold for the perturbed system. 

Similarly one applies an extended version of Theorem 11 

{ Tr \ ^ It 

by letting g ' (t ,x) = y in order to obtain a 
negative manifold for the perturbed system. 



CONCLUDING REMARKS 


The results of this dissertation may be considered 
as a parallel development of the works of several 
authors, [3], [12] , [15] and In particular they closely 
resemble Theorem 4,1, p, 330 In Coddington and Levinson 
[7] . These authors obtain results similar to Theorem 2 
for a system whose unperturbed state is linear, e . g. , 
system (4) in the absence of the 0-equation. Two aspects 
of the presentation in Coddington and Levinson are not 
given here. First, is is shown that solutions which do 
not start on the manifold eventually leave a neighborhood 
of the origin. Also it is shown that if the perturbation 
function is differentiable, as is the case here, then 
the manifold is likewise differentiable, i.e., the 
function which describes the manifold is differentiable. 
Although the manifold in the nonlinear case here has 
been shown to be Lipschitz, the problem of showing 
differentiability seems more difficult than in the 
linear case . 

Another aspect of the results in this dissertation 

which are not pointed out explicitly in the theorems is 

the similarity of behavior of solutions of the unper- 

+ + 

turbed and perturbed systems on S and S £ , respectively. 
Let (t , 9 (t ) ,x(t ) ,0 ) be in S + and (t ,if>(t ) , £ (t ) , w(t ) ) be 
in S*. Then j x (t )-£ (t ) | and |w(t)| approach zero as t 



becomes large. In the absence of the 0 -equation this 
may be interpreted as asymptotic equivalence of systems 
(8) and (11). This resembles results obtained by Marlin 
and Struble [19] who .give sufficient conditions for the 
existence of solutions for a perturbed nonlinear system 
which are asymptotic to solutions of the unperturbed 
state. Their hypotheses do not cover the case treated 
here, nor conversely. 
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ABSTRACT 


The purpose of this dissertation is to give a proof 
for the existence of an integral manifold for a system 
of perturbed nonlinear differential equations in a 
neighborhood of a critical point, periodic orbit, or 
periodic surface. Analogous studies of integral 
manifolds where the unperturbed system is linear have 
been done before. Functions describing the manifolds, 
in the linear case, are usually obtained as solutions 
of a certain improper integral equation formulated by 
the use of the classical variation of constants tech- 
nique. To prove that a solution to this integral 
equation exists, use is made of certain exponential 
bounds induced by the unperturbed linear system. Of 
interest here is the case in which the unperturbed state 
of the system has no linear part . Again one is led to 
consider solutions of an improper integral equation, 
obtained in this case from a generalization of the 
variation of constants formula due to V. M. Alekseev. 

In general, one can not expect exponential bounds 
analogous to those present in the linear case; however, 
by assuming certain smoothness and order type properties 
of the functions involved it is possible to demonstrate 
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the existence of a unique solution of the integral 
equation. Additional work is done to show that the 
solution of the integral equation gives rise to a 
solution of the differntial equation. 

It is also shown that on the manifold of the 
perturbed system solutions are asymptotic to solutions 
of the unperturbed system on the corresponding manifold. 



